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Abstract 


A developable surface can be formed by bending or rolling a planar surface without stretch- 
ing or tearing; in other words, it can be developed or unrolled isometrically onto a plane. 
Developable surfaces are widely used in manufacturing with materials that are not amenable 
to stretching. A ship hull design entirely composed of developable surfaces would greatly 
reduce production costs of that hull. 

This thesis describes a new, user-friendly method of designing developable surfaces with 
a B-spline representation. First, it expands the work of Aumann in designing developable 
B-spline strip surfaces whose directrices lie on parallel planes. The computer program 
developed can assist in designing surfaces from curves of any degree with any number of 
segments, and includes a test for regularity. Second, a new method is developed which 
permits the design of developable surfaces with general three-dimensional space curves as 
directrices. The computer program permits design of degree (3-1) B-spline developable 
surfaces with two patches using a minimization process. The basis is provided to extend 
this to surfaces with more patches and higher degrees. A test is provided to ensure that the 
minimization process results in a developable or nearly developable surface. This method 
is extended to include special interesting cases such as triangular degenerate patches and 
surfaces with a planar resulting directrix. 

This thesis also treats common differential geometry properties such as lines of curvature 
and geodesics that are useful in the design and manufacturing process. Lines of curvature 
are vital in the forming process, since the planar shape must be placed so the rollers are 
parallel to the lins:. o* sero curvature. As an inflection line greatly affects the forming 
process, a method is aescribed to determine the inflection line in advance; several properties 
of developable surfaces related to inflection lines are also described. 

Straight lines on a plane map to geodesics on a developable surface. This fact is of 
assistance in many aspects of the use of developable surfaces, including layout and quality 
control. A method is described to calculate geodesics on a developable surface as an initial 
value problem instead of the more complicated boundary value problem. 

Engineering examples are provided for each topic, and a small boat hull is designed 
using the methods described. In addition, recommendations for further research are made. 
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Chapter 1 


Introduction 


1.1 Background and Motivation 


A ruled surface is a curved surface which can be generated by the continuous motion of 
a straight line in space along a space curve called a directrix. This straight line is called 
a generator, or ruling, of the surface. Developable surfaces are a subset of ruled surfaces 
which have the same tangent plane at all points along the generator. A developable surface 
can be formed by bending or rolling a planar surface without stretching or tearing; in other 
words, it can be developed or unrolled isometrically onto a plane. Developable surfaces are 
also known as singly curved surfaces, since one of their principal curvatures is zero. 


Developable surfaces are widely used in manufacturing with materials that are not 
amenable to stretching. Applications include the formation of ship hulls, ducts, shoes, 
clothing and automobile parts such as upholstery and body panels [14]. 


In shipbuilding, developable surfaces are shaped using only rollers or presses. Heat 
treatment is then used only to remove distortion induced by welding or other means. Dou- 
bly curved surfaces, on the other hand, must be heat treated after rolling to induce the 
additional curvature. The heat treatment is normally done by hand, by a skilled artisan 
with years of training to achieve the correct amount of bending. This is an extremely time 
consuming, labor intensive and thus expensive process. 


According to Avondale/IHI Shipbuilding Technology Transfer data for a tanker, only 
15.1% of the curved plates in a ship hull are singly curved, while 65.8% of the plates are 
doubly curved, requiring roller and heat treatment processes [28]. An intensive effort to 
increase the amount of developable surfaces in the hulls of merchant ships at Burmeister 
& Wain Shipyard [34] has resulted in a reported 20% reduction in manhours required to 
produce a hull. Designing a ship entirely of singly curved, or developable, surfaces would 
reduce construction costs even more. 


Recently researchers in Computer Aided Geometric Design have been quite active in 
investigating the representation of developable surfaces in terms of Non-Uniform Rational 
B-Spline (NURBS) surfaces or a special case of NURBS called Bézier surfaces {1, 3, 14, 
29]. NURBS curves and surface patches are the most popular representation method in 
CAD/CAM due to their generality, excellent properties and incorporation in international 
standards such as IGES (Initial Graphics Exchange Specification) and STEP (Standard for 
the Exchange of Product Model Data). Thus, it would be beneficial to design developable 
surfaces using a B-spline representation. 


1] 


1.2 Research Objective 


Recently, efforts have been underway to revitalize commercial shipbuilding in the United 
States. Several professors at MIT are involved with developing cost saving methods in ship 
fabrication to ensure that shipbuilding in the US will be competitive in the world market. 
This research is in support of one of those efforts. 

The main goal of this research is to develop a user-friendly method of designing devel- 
opable surfaces with a B-spline representation. The effort is then extended to address some 
common differential geometry properties that will be useful in the design and manufacturing 
process. 

The ultimate goal is to provide a method to design a complete ship hull from developable 
surfaces and to generate cutting and bending information in a format that is user friendly 
for both the engineer and the worker. Although this thesis does not go that far, it takes a 
major step toward this goal. 


1.3 Thesis Organization 


In this thesis, Chapters 2 through 4 review basic differential geometry properties and intro- 
duce the concepts of developable surfaces. They also briefly review the Bézier and B-spline 
representations of curves and surfaces. The information introduced in these chapters will 
be used throughout the thesis. 

Chapter 5 describes a new user friendly method for design of developable surfaces in 
B-spline representation. Section 5.1 reviews the current literature on developable surfaces. 
Section 5.2 describes the design of strip surfaces that are constrained between two parallel 
planes. Section 5.3 describes the design of surfaces with directrices that are space curves 
and treats special cases such as triangular degenerate patches and surfaces with a planar 
resultant curve. Section 5.4 describes the unrolling, or development, of the surfaces into a 
plane. Section 5.5 gives some engineering examples. 

In Chapter 6, a geodesic on a developable surface is found as the solution to an initial 
value problem rather than a two point boundary value problem. The geodesics can be used 
in ship design for laying out butts and strakes, among other uses. In Chapter 7, the lines of 
curvature on a developable surface are analyzed and the line of inflection is defined. These 
concepts are required for the bending of steel plates into developable surfaces; the steel 
must enter the rolls in a direction parallel to the lines of zero curvature, and cannot be 
rolled past a line of inflection [33]. 

A small boat was designed using the methods described in this thesis. The boat is 
presented in Chapter 8 as a practical engineering example. 

The thesis concludes with Chapter 9 which includes recommendations for further re- 
search. 
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Chapter 2 


Review of Curves and Surfaces 


2.1 Introduction 


This chapter reviews the basic theory of curves and surfaces including such topics as the 
Serret-Frenet formulae, the first and second fundamental forms, curvatures and geodesic 
curves which will be used in derivations of many formulae in the later chapters. The 
information in this section follows the derivations in classical texts on differential geometry 
and geometric modeling such as those by Kreyszig [26], Struik [41] and Faux and Pratt [11]. 


2.2 Basic Theory of Curves 


Throughout this thesis, curves are represented in a parametric manner, r(u), as a function 
of one parameter, u, that lies within a closed interval uj <u < ug. The curve is a mapping 
from a one parameter interval to three-dimensional space 


r(u) =[z(u) y(u) 2(u)]’. 


The points on a curve are regular as long as at least one of the first derivatives is not equal 
to zero. In other words, a point is singular if dr/du = dy/du = dz/du = 0. The curve can 
be reparameterized if u can be expressed as u = f(uj) as long as du/du; # 0. 

Arc length is defined as the distance along a curve between parameter values ug and u, 
and can be represented as 


s(u) = | Vit+ P42 du= / Vr-r du. (2a) 
uo uo 


In this thesis, derivatives with respect to the arc length s will be represented by a prime 
and derivatives with respect to u will be represented by a dot. Derivatives of arc length s 
with respect to parameter u and vice versa are 





ds — 

oa bear — |r! 
Gs eaer er 

ed Wt 

ot oe 

U = ds |r| 
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tt du! r . r 


ds (&-#)?” 





The unit tangent vector t to a curve is defined by the unit vector that passes through 
two points in the curve, r(s) and r(s +h) (see Figure 2-1), as h approaches zero; i.e. 


BS 2) (see ee 








a es h aus! 
_ dr(u)du _ r(u) 
du ds |r(u)| 


This final representation clearly shows that t is a unit vector. The line passing through the 
point r(u,) in the direction of t(un) is the tangent to the curve at up. 





Figure 2-1: Tangent Vector 


All vectors normal to t(un) at up lie in a plane called the normal plane. Much as 
the tangent is defined by two points that approach one another on a curve, the osculating 
plane is defined by three points that approach one another. In other words, the osculating 
plane is the plane in which the curve lies at a point on the curve. If the curve is planar, 
the entire curve lies in a single osculating plane, and the osculating plane is constant for 
the entire curve. To mathematically define the osculating plane following the derivation in 
Kreyszig (26, p.31], take three points at parameter values u, u +h; and u+ ho. The plane 
can be defined by the two linearly independent vectors a; = r(u+h;) —r(u), «1 = 1,2 or 
linear combinations thereof. Let 


y (i) 2 r(u =f h;) i r(u) 
hj 
and 


2(v(2) — v(1)) 
hg — hy 


we 


14 


The Taylor series expansion of r(u + h;) is 


r(ut+h;) =r(u) + Air(u) + Meu) + 0(h?) (22) 


where o(h”) is a vector of Landau symbols o(h”) with the property that [41, p.3] 





Hay 
h>0 Ah” 
Using (2.2), we obtain 
=e r(u + h,) — r(u) 
hy 
_ r(u) + hit(u) + 5h7#(u) + o(h7) — r(u) 
—) = a ae 
2 
= r(u)+ Sera) + atin 
2, hy 


and 


2(v') _ vit) 
ho — hy 
2 
a = i, (#0 a -hai(u) + ho - Gc + haF(u) te nD 
Bay 2(h10(h3) — h2o(ht)) 
hyho(he—hi) 











As h, > 0, v") > r(u) and w > r(u). Therefore, the osculating plane is defined by r(u) 
and r(u). Note that r(w) must be at least twice continuously differentiable. 


osculating plane 


Figure 2-2: Osculating Plane 


The principal normal is defined as the intersection of the osculating and normal planes. 
Using an arc length parameterization, we find that differentiating t-t = 1 with respect to 
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s yields 


so t’ is perpendicular to t and lies in the normal plane. Since t = r’ = ru’, 
t! = ru" + i(u')*. 


Thus t’ lies in the plane of r and r, and therefore in the osculating plane. Since t’ lies in 
both the osculating and normal planes, a unit vector in the direction of t’ is the unit normal 
vector 
‘ 
n= —. 
jt"| 

The magnitude of t’ is called the curvature and is represented as k. The curvature vector 
k is represented as 


k= t = An. (253) 


A physical sense of the curvature can be found from the simple example of a line, which 
has a constant tangent vector, so t’ is zero. Therefore, « = |t’| = 0, and a line has zero 
curvature. Parameterized by u, k and « can be represented as [41, p.16] 


Lee a eye “a =! tr 
2 Coan eet), (eon ae Li <7) 
ae oe | | (r-£r)7 (r-£)3 


where the Lagrange identity (a x b)- (a x b) = (a-a)(b- b) — (a- b)? is used. 
The normal to the osculating plane is the binormal, which is defined as 
b=t Xn. 


Since these three vectors are unit vectors that are mutually orthogonal, they satisfy the 
relations 


Differentiating b- b = 1 and b- t = 0 yields 
b’-b=0 
and 
b'-t=—b-t’ = —b- «n= —x«(b-n) =0. 
Therefore, b’ is orthogonal to both b and t and must then be parallel to n. Let 


b’ = —7n 
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where 7 is the torsion of the curve. A physical sense of torsion can be found from the fact 
that a curve with zero torsion and nonzero curvature is a planar curve. 


Finally, an equation for the derivative of the normal vector corresponding to the rate of 
change of the tangent plane (26, p.41] can be determined. Differentiating n-n = 1 yields 
n-n’ =0. Therefore, n’ is orthogonal to n and must satisfy the relationship 


n’ = at + Bb. (2.4) 

Multiplying equation (2.4) by t and b respectively yields 

Y= t ands —n bp. 
Differentiating n-t = 0 yields 

a=n’-t=-—-n-t’=—-n-«Kn=-—k«. 
Similarly, differentiating n- b = 0 yields 
B=n'-b=-—n-b’ = —n-(-7n) =7. 

Thus, 


n’ = —Kt + 7b. 


The representations of the derivatives of the three defining vectors for a curve are termed 
the Serret-Frenet formulae {41, p.18] 


t = «nN 
b = -7rn 
no = —Kt+7Tb, 


which describe all aspects of a space curve, as shown in Figure 2-3. 


Osculating Plane 





n 
Normal! Plane 5 os 


Tangent Plane 





Figure 2-3: Space Curve 
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2.3 Basic Theory of Surfaces 


A general parametric surface can be defined as a vector-valued mapping from a two- 
dimensional parametric uv-space to three-dimensional space 


R(u,v) = (z(u,v) y(uv) z(u,v)]", 


where the two parameters u and v lie in the closed intervals u,; < u < ug and v; < v < v2. 
The surface is regular at all points where the matrix 


Or Oy Oz 

=o OQ OQ 
M = oe (2.5) 

Ov Ov Ov 





has rank 2, or equivalently R, x R, 4 0 [41, p.56]. Partial derivatives of a surface will be 
represented as 
OR OR 


Ba = and = a 


A curve on the surface is represented as r(u, v) with u = u(t) and v = v(t) or equivalently 
R(u(t), v(t)) = r(t). This occurs when the rank of the matrix (2.5) is 1 at every point. 
When we keep u constant in R(u,v) by setting u = un, we have an equation for a curve 
that depends only on v, and is thus called an isoparameter curve. Similarly, v = vp is also 
an isoparameter curve. 


The vector 


dR _ 
dt 


du du 
a ae Lh 2.6 
Rua + Roa, (2.6) 
is tangent to the curve when wu and v are functions of ¢t. This is of course also a tangent to 
the surface. The tangent plane to the surface at any point is determined by the two vectors 
R, and R,. The normal to the surface is orthogonal to this plane. A unit vector in the 


direction of the surface normal is given by 


R, x R, 
IR, x Ry| 


Note that this definition of the normal vector is undefined when R, x R, = 0, supporting 
the assertion that a point on the surface is singular if R, x R, = 0. This may be due to 
the shape of the surface or the choice of parameterization. In Section 5.3.5, an alternative 
definition of N is explored for degenerate points. 


Nie (2am 


Using the definition of arc length (2.1), the distance between two points on a curve on 
the surface is found by integrating 


ds = Vdzrdz + dydy + dzdz = VdR- dR. (2.8) 
Combining (2.6) and (2.8) yields the first fundamental form of the surface 


I = ds? = (R,du + R,dv) - (R,du+ R,dv) = Edu? + 2Fdudv + Gdv? 
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where 
E=R,-R,, F=R,-Ry, G=R,-Ry. 
Note that 


EG 


(R,,-R,)(R,-R,) — (R,-R,)? 
SE eet) EtG) 
(Ry, x R,)|?. (2.9) 


Thus, EG — F? is always positive if the surface is regular. 


The first fundamental form enables us to determine arc lengths, angles and areas on the 
surface. In order to determine curvatures, the second fundamental form is required. First, 
the curvature vector of the curve on the surface, found from equation (2.3), is split into its 
normal and tangential components 


le =F = en + hey = aN + ky 


where k,, is the normal curvature vector and Kk, is the normal curvature of the curve on 
the surface at that point. The tangential portion of the curvature vector, Kg, is called the 
geodesic curvature vector. It will be discussed in more detail later. 


If the angle between the surface normal and the normal to the curve is designated + as 
shown in Figure 2-4, then [26, p.118] 


n- N = cos(7)|n||N| = cos(7) 
from the definition of a dot product. Since t’ = R” =«-n, 


Kcos(y) = «-n-N 

R” : N 

(R,du + R,dv)'-N 

(R.,,,dudu + 2R,,dudu + R,y,dvdv + R,d?u+R,d*v)-N. (2.10) 


Since R,-N = R,-N = 0, equation (2.10) reduces to the second fundamental form, 
represented by 


II = (Ryydudu + 2R,,dudu + Ryy,dudv) -N 


= Ldudu+2Mdudv + Ndvudvu 
where 
L Riu: N (Zea 
i ites Ny) and (Dal2) 
N = R,,-N. (2.13) 
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Also note that since R,-N =0, 


a = Ri,:-N+R,-N, =0, 


and therefore 
Rye N = —-R,-Ny. 


Similarly, Ry, -N = —R,-N, = —R,-N, and R,,-N = —R,-N,. Therefore, 


iL = =RaaNe (2.14) 
M. = {RNR (2.15) 
= -5(Ry-N,+R,-N,) and 
N 





Figure 2-4: Curvature Vectors 


To find the normal curvature, we begin with the fact that the surface normal is orthog- 
onal to the tangent to the curve. Differentiating t -N = 0 with respect to s yields 


dt dN 
Thus we have 
dt dN dR dN 
ar ds deems 
Therefore, 
dR dN (R,du+R,dv) - (N,du + N,dv) JOE 
Kk, = —-k-N= —.: — = eo, = - 


ds_ ds ds? I 
Ldudu + 2Mdudv + Ndvdv 


~ Edudu + 2Fdudv + Gdudu 
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_ _L+2MA+ NX (2.17) 
—  B+2FX +62 


where A = du/du. The negative sign in equation (2.17) ensures that if kK, is positive, the 
center of curvature lies opposite to the direction of the surface normal. 


Any point on the surface has many curves passing through it, each of which has an 
associated k,. The maximum and minimum values of ky, at a point are called the principal 
curvatures of the surface at that point and are designated «, and K2. The maximum and 
minimum values of «, occur when Se =O iene ln) 


on = (E+ 2FX+G)*)(2M + 2Nd) —(L+2MA+ N*)(2F +2G) =0. (2.18) 


Using (2.18) in (2.17) yields 


«. - EtMA__M+NA 
ree nee OG 


given 


(EB + 2F) + Gd’) 
(L+2Md + Nd?) 


(B+ FA)+A(F + GA) and 
(L+ MA) +A(M+N)). 


Therefore the extreme values of kK, satisfy the two simultaneous equations 


(L+knB)du+(M+n,F)dv = 0 and (2.19) 
(M+ky,F)du+(N+kn,G)dv = 0. (2.20) 


These equations can be simultaneously satisfied if and only if 


L+tknb M+e«,F 


MoU ere Glee 


The discriminant of this quadratic equation in «,, is greater than or equal to zero. Therefore 
the equation has either two distinct roots Kmez and Kmin, the maximum and minimum 
principal curvatures, or double roots «,, the normal curvature at an umbilical point. The 
corresponding directions define directions in the uv-plane. The corresponding directions 
in the tangent plane are called principal directions of curvature and are always orthogonal 
except at umbilical points, where Kmax and Kmjn are identical. In the special case where the 
identical principal curvatures vanish, the surface becomes locally flat. The two roots are 
given by 


Kmxe = H+W7H?-—K and (2.21) 
Kmin = H-WJH?-K (2.22) 


where K is the Gaussian curvature and H is the mean curvature defined by 


LN — M? 
eC = ee | 
EG. FE (2723) 
2hPM—EN—GL 
- “3EG-F) ii 


a 


Using (2.21) and (2.22) in (2.23) and (2.24), it can easily be shown that 


K = Kmoarkmin 
1 
Hq = 9 (Kmar ats Forni 


Note that at an umbilical point, H? = K and at a flat point, K = H = 0. 


The normal curvature can be represented in terms of the principal curvatures K, and K2 
-as 


Kn = 1 coS?(a) + Ke sin*(a) (2.25) 


where a is the angle between the direction du/du and the direction dv = 0. This is known 
as Euler’s theorem [26]. 


Recall that the curvature vector of a curve on the surface can be split into its normal 
and tangential components 


k =k, + ky, 


where the tangential component, k,, is known as the geodesic curvature vector. Let u be 
a unit vector in the tangent plane perpendicular to t such that t, u and N form a right- 
handed coordinate system in that order as shown in Figure 2-4. The geodesic curvature Kg 
can then be found from 


k, = KgU. 
Additionally, since the magnitude of «, can be represented as | cos y|, the magnitude of 
Kg can be represented as |K sin 7]. 


Geodesic curvature can be represented using only the first fundamental form. Since 
u-N=0O0andu-u=l1, 


dt 
eet 7 =u t = (Nxt) 6 —(tt N). (2.26) 


The unit tangent vector of the curve on the surface is given by 
t=R,u’ + Rv’; 
thus 
t’ = Ry (u’)? + 2R,yu'v’ + Ryy(v')? + Ryu" + Rv". 
Therefore equation (2.26) can be rewritten as 


Ko = (( xian 
(Ryu + Ryv’) x (Ruu(u’)? + 2Ryyu'v’ + Ryy(v’)? + Ryu” + Ryv")]-N 
(Ra, x Ruw)(u’)? + (2Ry x Ruy + Re x Ruy) (u’)? 0’ 
+(Ry x Ryy + 2Ry x Ryy)u'(v’)? 
+(Ry x Ryy)(v')? + (Ry x Ry)(u'v” — uv’) .N. O27 
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Each of the coefficients of the combinations of u and v derivatives in (2.27) can be written 
in terms of the coefficients of the first fundamental form and their derivatives. For example, 


(Re foo S§GoFD 
VEG — F? 
2(E (Raw > R,) — (Ra: Run) F)V EG — F? 


2(EG — F2) 
_ (EF, — EE, — FE,)VEG — F? 
7 2(EG — F?) 


= I? VEG - F? 


where I°?, is one of the Christoffel symbols defined by the coefficients of the first fundamental 
form E, F and G and their derivatives, E,, Fy, Gu, Ey, Fy and G, 


GE, —2FF,+ FE, 


2 ee 
eae 


2(EG — F?) 2(EG — F?) 
ri — Gh FG. pe, . BGu-FE, 
on en G ere a SO(n Gree) 
a HEG-—F2) ’ 42° s« O(EG = F?) 
Similarly reducing all the other coefficients in (2.27) yields 
a (Pt, (u')° ie (2Pfo a Di1)(u’)?0" na (a am 20 j2)u'(v')? - ['39(v')? 


tu'v” —ul"v'|VEG — F?. (2.28) 


A geodesic is defined as a curve with zero geodesic curvature [41]. Straight lines on a 
surface are geodesics since the curvature vector k vanishes. For geodesics that are curved, 
the curvature vector coincides with the surface normal vector and, since the curvature vector 
lies in the osculating plane, the osculating and tangent planes are normal. The equation for 
the geodesic can be obtained by setting «, equal to zero in equation (2.28), which yields 


uy! —ulo" = Thy (w)? + (2%, — Ph) (u!)?0" + (C9 - Vha)u(v')? — Th (o")? 


assuming that the geodesic is everywhere regular so EG — F? is always positive. An al- 
ternative representation for a geodesic is given by the set of coupled second order ordinary 
differential equations [41] 


d7u du \ 7 du dv dv \? 

ais (C2 eo ce »(S) 2 2.2 
oe h($) Sie a2 eee 2n\ es oa oe») 
d2y du \ 2 du dv du \ 2 

me ope ee (=) Z . 
ds? * i (S) x M12 Gs ann ds ° et 


where the two equations are related by the condition ds? = Edu? + 2Fdudv + Gdv?. 
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Chapter 3 


Curve and Surface Representation 


3.1 Introduction 


This chapter reviews the representation of curves and surfaces in Bézier and B-spline forms 
and treats the special properties associated with each. In addition, algorithms for knot 
insertion and removal, curve evaluation and splitting are summarized. The descriptions are 
based on texts by Hoschek and Lasser [22], Piegl and Tiller [36] and Yamaguchi [44]. 


3.2 Bézier Curves and Surfaces 


Bernstein Polynomials 


The Bernstein polynomials are defined as 


_ 
nett a be yet Owen (3.1) 


7 ilin Dl 


and have several properties of interest. The property of positivity states that each polynomial 
factor is non-negative such that 


Bin(u)>0, O<us<il 


for alli and n. The partition of unity property states that the Bernstein polynomials sum 
to 1 for allO <u<1,or 


nm 


> Bn) ae 
7=0 


The derivative of a Bernstein polynomial is 


aPinl) = RB ernst) a ee 


where B_j n-1 = Ban-1 = 0. The linear precision property of a Bernstein polynomial 
eal? 
C= S- 7 Bin(u) 
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allows the monomial u to be expressed as the weighted sum of Bernstein polynomials with 
coefficients evenly spaced in the interval [0,1]. 


Bézier Curves 


A Bézier curve is a spline curve that uses the Bernstein polynomials as a basis. A Bézier 
curve of degree n (order n + 1) is represented by 


Oa <= 1. 


a 
= 

! 
M4 
s 
wD 
3 
= 


t=) 


The coefficients, b;, are the control points that determine the shape of the curve. Lines 
drawn between consecutive control points of the curve form the control polygon. A cubic 
Bézier curve is shown in Figure 3-1. Bézier curves have the following properties: 


e The first and last control points are the endpoints of the curve. In other words, 
bo = r(0) and b, = r(1). 


e The curve is tangent to the control polygon at the endpoints. 


e The conver hull property states that the entire curve is contained within the convex 
hull of the control points. 


e The first derivative of a Bézier curve is represented by 





r(u) = =n) (bi41—b)Bin-i(u), OSu<l. 





Control Polygon oe is ys 


Figure 3-1: A Cubic Bézier Curve with Control Polygon 


A Bézier curve can be evaluated at a specific parameter value wp and the curve can be 
split at that value using the deCasteljau algorithm, where [22, p.125] 


bi (uo) = (1 — up) be) + ub; ', &£=1,2,...,n, i=k,...,n 
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is applied recursively to obtain the new control points. The algorithm is illustrated in 
Figure 3-2, and has the following properties: 


e The values b? are the original control points of the curve. 
e The value of the curve at parameter value ug is b7. 


e The curve can be represented as two curves, with control points (bj, b7,..., b?) and 
ipee br ',..., b2). 


e Since this is merely a change in the basis representation, the curve itself remains 
parametrically and geometrically unchanged. 





Figure 3-2: The deCasteljau Algorithm 


Bézier Surfaces 


A tensor product surface is formed by moving a curve through space while allowing defor- 
mations in that curve. This can be thought of as allowing each control point b; to sweep a 
curve in space. If this surface is represented using Bernstein polynomials, a Bézier surface 
is formed, with the following formula: 


m mn 
R(u, v) i= 3 S- bi; Bim(u)Bjn(v), ess le 
7=0 j=0 


Here, the set of lines drawn between consecutive control points bj; is referred to as the 
control net. An example of a bi-quadratic Bézier surface with its control net can be seen in 
Figure 3-3. The following conditions apply: 


e The boundary isoparametric curves (u = 0, u = 1, v = 0 and v = 1) have the same 
control points as the corresponding boundary points on the net. 


e The corners of the surface coincide with the corner points of the net, and the deriva- 
tives in the u and v directions at the corners are tangent to the net. 


e Bézier surfaces exhibit the convex hull property. 
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Figure 3-3: A Bézier Surface with Control Net 


Continuity Conditions 


Bézier curves can represent complex curves by increasing the degree and thus the number 
of control points. Alternatively, complex curves can be represented using composite curves, 
which can be formed by joining several Bézier curves end to end. If this method is adopted, 
the continuity between consecutive curves must be addressed. 

One set of continuity conditions are the geometric continuity conditions, designated by 
the letter G with an integer exponent. Position continuity, or G° continuity, requires the 
endpoints of the two curves to coincide, 


r)(1) = r'2)(0). 


The superscripts denote the first and second curves. Tangent continuity, or G! continuity, 
requires the tangents of the curves to be in the same direction, 


r(1) =ait 
r'2)(Q) = aot 


where t is the common unit tangent vector. G! continuity is important in minimizing stress 
concentrations and preventing flow separation. 

Curvature continuity, or G* continuity, requires the center of curvature to move contin- 
uously past the connection point, 


z2)(0) = (22) 2) + pr")(1). 


G* continuity is important for aesthetic reasons and for preventing fluid flow separation. 

More stringent continuity conditions are the parametric continuity conditions, where 
C* continuity requires the kth derivative of each curve to be equal at the joining point. In 
other words, 


d¥re(1) — d*r'2)(0) 
duk ——s dusk 
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The C! and C? continuity conditions for consecutive segments of a composite degree n 
Bézier curve can be stated as [44, p. 215] 


eae (Dni = Br=1) = h; (apie a bni) ‘ and fo. 2) 
h; De 
bni-1 + a2 (been ~— bni-2) = Dniti + ——(Dnigi — Dni+2) (ais) 
h; Ais) 


where, for the ith Bézier curve segment with u values between knot values ¢; and t+ 
inclusively, h; = tj+, — tj. See Figure 3-4. 


by, +3 





Figure 3-4: Continuity Conditions 


3.3  B-Spline Curves and Surfaces 


The Bézier representation has two main disadvantages. First, the number of control points 
is directly related to the degree. Therefore, to increase the complexity of the shape of the 
curve by adding control points requires increasing the degree of the curve or satisfying the 
continuity conditions between consecutive segments of a composite curve. Second, changing 
any one control point affects the entire curve or surface, making design of specific sections 
very difficult. These disadvantages are remedied with the introduction of the B-spline 
representation. 

A B-spline (basis-spline) curve is a piecewise polynomial curve meaning it is made up of 
polynomial curve segments. It is a spline curve with a different basis than the Bézier curves, 
although a Bézier curve is a special case of a B-spline curve. The B-spline representation 
has the advantage that the control points affect the curve only locally, and a B-spline of 
order k may have as many control points as required to describe the curve. 


A B-spline curve is a series of polynomial segments of degree k ~— 1 joined together at 
knots t,. A B-spline of order k (degree k ~ 1) with n+ 1 control points has the equation 


r(u) = PiNix(u), eels a € [f.-1, tn41| 
Th 
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where t; is a member of the knot vector with n+k+1 elements 
T= (to, bi, -++y bei, bk, tk+1) seeyln—-1,tn, En+1; oeehsiy tn+k)s Comet es tn+k: 
The basis functions, N;,(u), are given by: 


1 for in <u< Git 
Q otherwise 


Nj1(u) = 


for k = 1, and 


u—t; Cea 


Ne 
iw(U) bitk — ti41 


Niz—-i(u) + eine) 


fork > 1 andi =0,1,...,n. These equations have the following properties [22, p. 168]: 
oe (tt) 0, tote ee aire 
e Nix(u) = 0, for to <u < tj, and ti4, Su < tne 
oe Nor) 4, for w eteninta and 
e Nx (u) has continuity C*~? at each simple knot ty. 


The first three properties show positivity and partition of unity for the B-Spline basis 
functions. 


“Pr 





Figure 3-5: A Cubic B-spline Curve Segment with Control Polygon 


Local control means that a single segment of the B-spline curve is controlled only by the 
nearest n points, and that any control point affects the nearest n spans. In other words, 
changing p,; affects the curve in the parameter range t; < u < t;,, and the curve at a point 
u where t, < u < t,4) is determined completely by the control points p,_(,4~1),---, Pr- 

The convex hull property for B-splines applies locally, so that a span lies within the 
convex hull of the control points that affect it. This provides a tighter convex hull property 
than that of a Bézier curve, as can be seen in Figure 3-6. 
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Figure 3-6: Convex Hull Property for a Cubic B-Spline Curve 


Increasing the multiplicity of a knot reduces the continuity of the curve at that knot. 
Specifically, the curve is (k —p—1) times continuously differentiable at a knot with multiplic- 
ity p (< k), and thus has C(#-?-)) continuity. Therefore, the control polygon will coincide 
with the curve at a knot of multiplicity k — 1, and a knot with multiplicity k indicates C~! 
continuity, or a discontinuous curve. 


As you can see in Figure 3-5, the endpoints of the curve may not coincide with the 
control polygon. However, repeating the knots at the end k times will force the endpoints 
to coincide with the control polygon. Thus a knot vector described by 


een ths bky bkt1)-++)'ln—-1, En, oscil +> vaaleemet ke) 
a Ee 


k equal knots n—k+1 internal knots k equal knots 


Pal@iave tg = t) =... = tp_) and tas) = th42 = :.. = typo, and hence the first and last 
control points of the curve coincide with the endpoints of the curve. 


From this discussion, it can be seen that a Bézier curve of order k (degree k —1) is a 
B-spline curve with no internal knots and the end knots repeated k times. The knot vector 
is thus 


T= (to; ei e b eins Caen) 
——_— Oe” 


k equal knots &k equal knots 


wheren+k+1=2k orn=k-1. 
A B-spline surface is a tensor product surface using a B-spline basis. This is represented 
mathematically as 


m n 


R(u,v) = >>> pi Nix(u)Nj,(v). 


1—0 7=0 
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Knot Insertion 


A knot can be inserted into a B-spline curve without changing the geometry of the curve. 
The new curve is identical to the old one, with a new basis where 


n . n+1 - 
S- piNie(u) becomes De PiNit (u) 
1=0 1=0 
overly — lip, ti... eee Guer i = (to, f),..., t,t, fee 


when a new knot ¢ is inserted between knots t, and ti4,. The deBoor algorithm is used to 
insert the knot, and is represented as 


p]=(1—a/)pip+ajpp'  izl—k+2 ee 
where 
1 i<l—-k+1 
a, = ¢ 0 l+1<i 


pe=k 1-k+2<iK<! 

[+k~—1 7% 

To evaluate a B-spline curve at a specific parametric value or to split the curve at that 
value, the deBoor algorithm is used to repeatedly insert the same knot until the control 
polygon coincides with the curve at the parametric value uo, or in other words, until the 
multiplicity of the knot at w= k—-—1. Then, 


P; = pi 
Dar eeGe 


A B-spline curve is C™ continuous in the interior of a span. Inserting a knot does not 
change the curve, so it does not change the continuity. However, if any of the control points 
are moved after knot insertion, the continuity at the knot will become C*~?—!, where p is 
the multiplicity of the knot. Figure 3-7 illustrates a single insertion of a knot at parameter 
value ug, resulting in a knot with multiplicity of one. 


ae ederim een a5 





p3=p4 


Figure 3-7: The deBoor Algorithm 
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The B-spline curve can be subdivided into Bézier segments by knot insertion at each 
internal knot until the multiplicity of each internal knot is the degree k. The above algorithm 
is also known as Boehm’s algorithm [4, 5]. 


Knot Removal 


Knot removal is the reverse process of knot insertion. This thesis uses the knot removal 
algorithm developed by Tiller [42], in which a knot is removed if and only if it leaves the 
curve or surface geometrically and parametrically unchanged. 

To demonstrate the process, this example uses a cubic B-spline curve r(u) given by 
control points (p3,...,pg) and knot vector (to,...t19) where tp =... = t3 = 0, tg = ts = 
ts = 1 and t7 =... = tig = 2 as shown in Figure 3-8. As the basis functions only guarantee 
C® continuity at u = 1, the first derivative may or may not be continuous there. Using the 
C! continuity condition (3.2), the first derivative will be continuous if and only if 


(t7 — t4)(P3 — po) = (t6 — t3)(p4 — P3)- 
Since t4 = tg = U, 


u—t3 t7 —uU 
3 pi + 2 
t7 — t3 








t7 — t3 


Since p} = p} and p} = p}, 





u—t 
p? =03p3+(l—a3)p, a3= : (3.5) 


ty —t3 


Note that (3.5) is the same as the deBoor algorithm equation (3.4) for inserting a new knot 
at u = 1, although the points and knots are numbered differently in this example. 

A similar reasoning yields that fact that a knot u = 1 can be removed a second time if 
and only if the second derivative is continuous, yielding 


P> = a2p5 + (1 — a2)pj 
P3 = a3p3 + (1 — a3)p5 
a oe (3.6) 


a; = ———————_ 
bitp+2 — ty 


and a knot wu = 1 can be removed a third time if and only if the third derivative is continuous 
such that 


Pj = ap} + (1 —a1)p9 


p35 = a2p3 + (1 — a2)p} 
p3 = 03p3 + (1 —a3)p3 

ee, ale (3.7) 
titpt3 — ti 


Note that there are no unknowns in equation (3.5), one unknown, p, in equation (3.6) and 
two unknowns, p? and p3, in equation (3.7). 

For the knot removal process, first the right hand side of equation (3.5) is computed and 
compared to p3. If they are equal within a given tolerance, the knot and p3 are removed. 
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If the first knot removal is successful, equations (3.6) are solved for p3 and compared: 


2 _ Py—(l—a2)pj 
bor =e 
a2 
p? = P3 — 23p3 
2 1—a3 


If the two values for p? are the same, then the knot and control points p} and p} are 
removed and control point p3 is inserted. 

If the second knot removal is successful, the third step is to solve the first and third 
equations of (3.7) for 


3 _ Pi-—(l—a)p% 
0 —————— 
Q) 
ps = P3 — 3p} 
g l1—ag3 


The two values are then substituted into the second equation of (3.7). If the result is within 

tolerance of p3, then the knot is removed and control points p{, p3 and p$ are replaced by 
3 3 

p) and pp. 





Po =Po=P6=P5 Pg =Ps5=P4=P3 


lost ie os! 5 tgats ate t7atg atg =to 


Figure 3-8: Knot Removal 


This can be generalized to apply to any number of removals of any particular knot. For 
the nth removal, there will be a system of n equations with n — 1 unknowns. If n is even, 
two values of the final unknown control point will be calculated and compared. If they 
are within tolerance, the knot removal is successful. If n is odd, all new control points are 
computed and the final two are substituted into the middle equation. If the result is within 
the tolerance, the knot removal is successful. If the nth removal is successful, n control 
points will be replaced by n — 1 control points. 

Knot removal from a surface is performed on the m+ 1 rows or n+ 1 columns of control 
points. However, the knot removal is successful only if the knot can be successfully removed 
from each row or column. Therefore, the result must be checked for each row or column 
before any control points are removed. 
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NURBS 


The non-uniform rational B-spline (NURBS) curve is represented in a rational form 


io WibiNi«(u) 
Sy WN Ele) 


where w; > 0 is a weighting factor and N;,(u) are defined over a non-uniform knot vector. 
A NURBS surface patch can be represented as 


Ru, v) = Liao Lazo Wis dis Nie (w)Njulv) 


where w;; > 0 is a weighting factor. A uniform knot vector is one where the knots are 
evenly spaced. If the knots are not evenly spaced, the knot vector is non-uniform. Thus, 
any knot vector with repeated knots is non-uniform. The NURBS representation of curves 
and surfaces allows the exact representation of figures such as circles, conics, quadratics, 
and surfaces of revolution with rational profiles such as toril. If w; = 1 for all i or w;; = 1 
for all 2,7, this representation is the same as the integral non-uniform B-spline described 
previously. In this thesis, all the B-splines used are integral, but the development could be 
extended to include NURBS. 


leg) 
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Chapter 4. 


Review of Differential Geometry 
Properties of Developable Surfaces 


4.1 Introduction 


This chapter discusses the properties of ruled and developable surfaces. The representations 
and properties discussed herein are taken from classic texts by Struik [41], Kreyszig [26], 
Faux and Pratt [11] and Spivak [40]. 


4.2 Ruled and Developable Surfaces 


Ruled Surfaces 


A ruled surface is a surface generated by the motion of a straight line (a generator or ruling) 
through space [41]. A curve that passes through all the rulings of the surface is called a 
directriz. Any point on the surface can be expressed as 


R(u,v) = a(u) + vB(u) (4.1) 


where a@(u) is a directrix or base curve of the ruled surface and ((u) is a unit vector which 
gives the direction of the ruling at each point on the directrix. Alternatively, the surface 
can be represented as a ruling joining corresponding points on two space curves. This is 
represented by 


R(u,v) = (1—v)ra(u)+urg(u) O<u,v<l (4.2) 


where r4(u) and rg(u) are directrices, as shown in Figure 4-1. The two representations are 
identical if 


a(u)=ra(u) and G(u) =rg(u) —ra(u). (4.3) 


Throughout the thesis it is assumed that the constant u isoparametric lines correspond to 
the generators of the developable surface or, in other words, the straight line rulings are in 
the v direction. 
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directnx, ra(u) 


generator, ruling 





directrix, fg (u) 


bee 


Figure 4-1: A Ruled Surface 


Developable Surfaces - A Subset of Ruled Surfaces 


Developable surfaces are a subset of ruled surfaces that have a constant tangent plane along 
each generator. Since surface normals are orthogonal to the tangent plane and the tangent 
plane along a generator is constant, all normal vectors along a generator are parallel. This 
is shown in Figure 4-2. 


ee surface normals 
tangent plane. | 


generato 


Figure 4-2: A Developable Surface with Tangent Plane along a Ruling 


Given two points on a single generator that are selected at parameter values (ug, v;) 
and (ug, v2), the surface is developable if the tangent planes at these points coincide. The 
tangent planes are defined by the vectors R, and R,. Using the representation of a ruled 
surface shown in (4.2), we obtain 


dr 4(uo) Hee dr (uo) 





R, (uo, v1) = (1 — v4) a a (4.4) 
R,(uo, v2) = (1 — ug) Talo) a i (4.5) 
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R, (uo, v1) = Ry(uo, v2) = re(uo) — ra(uo). (4.6) 


If the tangent planes at the two points coincide, the vectors R,(uo,v1), Ru(uo,v2) and 
rp(uo) —ra(uo) must be coplanar. Therefore, the triple product must equal zero: 


IRu(uo,vi) Ru(uo,v2) (re(uo) — ra(uo))| = 0. (4.7) 
Let R,(uo, v1) = Ry and R, (uo, v2) = Ry2. Then 
Rui x Ry: = (RY, RZ, — R?,R¥,]é — [R7,R2, — R2,R2,]j + (RZ,R¥, — RY, RZ) (4.8) 
Expanding the first term of (4.8) by employing (4.4) and (4.5) yields 
Be fo Ri Rip an v1) 0% + uy ite ae?) eee verp] 
—((1 = v1 )r4 + vrgi|[(1 — ai + vers] 
= (l—w)(1— v2)r4r%y + vi (1 — v2) rary + ve(1 — 01) Arg + vivergrg — 
(1 = v;)(1 = v2)r4r%, + v1 (1 - vo)rth + vo(1 _ UP GAL + V1 Ver Bry 
= (v1, — VV. — V2 + UL V2)FArR + (ve — vive — v1. + v1 v2)F4 rR 


= (v2 —1)(T4rB — FArp). 


Similarly, 

Ri Rig — Ri Rig = (ve — v1) (rare — rare) 
and 

RE Rip — RY Rip = (v2 — v1) (EE — PED) 
Therefore, 


Rai X Rug = (v2 —v1)(Fa X rp) 


and thus (4.7) can be reduced to 
(rp(u) — ra(u)) x ra(u)-ra(u) =0 (4.9) 


for any u as long as v; # vg. Conversely, any surface that satisfies (4.9) must satisfy (4.7) 
and, therefore, the tangent planes must be the same and the surface must be developable. 
Thus, a surface is developable if and only if (4.9) is satisfied. Substituting (4.3) into (4.9) 
and using the fact that r4 x r4 = 0 yields the equivalent condition 


ax fp-B=0. (4.10) 


Developable Surfaces - Envelopes of Families of Planes 


An alternative representation of a developable surface is as an envelope of a family of planes. 
The concept of an envelope will be described here for a family of curves, then extended to 
planes. First, we represent a curve as r(u,qa) where wu is the curve parameter and a is the 
family parameter. In other words, each curve in the family is r(u,a,). There may be a 
curve re(qa) that is tangent to each curve in the family, as shown in Figure 4-3. This curve 
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would be the envelope of that family of curves [11]. 


rs r(u,a) 


rs r(u,@) 





r= r(u,a,) 


Figure 4-3: An Envelope of a Family of Curves 


Similarly, a surface can be the envelope of a family of surfaces. Let a one-parameter 
family of surfaces be represented by G(x,a@), where a is the family parameter, and assume 
that the surfaces are twice continuously differentiable. Also assume that consecutive surfaces 
intersect one another. The intersection is 


G(x,a)=0 G(x,at+h)=0 
where |h| is sufficiently small. This may also be represented by 


G(x,a +h) — G(x, a) 
h 


The limit of (4.11) as h tends to zero is OG(x,a)/Oa = 0. The set of points represented by 


= 0. (4.11) 


OG(x, a) 
Oa 


is the characteristic of the surface. If the characteristics of all the possible surfaces in the 
family form a surface, that surface is called the envelope. If the envelope exists, then at 
every point of the characteristic of a surface in the family, the tangent planes of the surface 
and of the envelope coincide. 


Gio a)— 0) =) 


A characteristic point is determined by the intersection of the characteristic with another 
surface in the family. These points of intersection are determined by 


OG (x, a) 


G(x,a) =0 Fo 


=) (GleS esa is) = 0) 


where |k| is sufficiently small. Using the Taylor series expansion 


OG (x, @) A k? G(x, a + 0k) 


G(x,a+k)=G(x,a)+k Aa aa. 


OraO= 1, 
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the third equation can be rewritten as 


0*G(x, a + 6k) 
Oa? 


Then, taking the limit as k approaches zero, the characteristic points of the surface are 


= 0. 


OG (x, @) 0°G(x, a) 

G(x,a) = 0 a= mee ge 
If the characteristic points exist and if they form a curve as a@ varies, the curve lies on the 
envelope of the family and is called the edge of regression of the envelope. If the envelope 
and the edge of regression exist, then, at every point of the edge of regression, the tangents 
of the edge of regression and the corresponding characteristics coincide. The characteristics 
then form two sheets corresponding to the tangents of the edge of regression in the positive 
and negative directions. These sheets form a cusp at the edge of regression, as shown in 
Figure 4-4. 





Second Sheet 
Characteristic 


Figure 4-4: Edge of Regression 


Let us now consider a family of planes and their derivatives which can be represented 


by 


G(x,a) = a(a)-x+c(a) = 
OG Oa Oc 
da da" Ga 

0G dra a7 


AE ay 


where a(q) is the unit normal vector to the corresponding plane. When all planes are paral- 
lel, no intersection exists and no characteristic curves are formed. When all planes intersect 
in a single line forming a pencil, the envelope is a line instead of a surface. Therefore, these 
two cases are excluded [26]. When the envelope exists, it is formed of the characteristics of 
the family of planes. Since planes intersect in straight lines, the envelope is a ruled surface. 
Furthermore, the tangent plane along a characteristic is constant since this is an envelope 
of planes, and hence the surface is developable. 


4] 


When the vectors a(q) are all parallel to a single plane P, then all the planes in the 
family are parallel to a vector normal to the plane P. Thus, the envelope of the family of 
planes forms a cylinder (except when they form a pencil which has been excluded). When 
all the characteristic points coincide, all the characteristics intersect in a single point and 
the surface is a cone. When the characteristic points form a curve, the curve is the edge 
of regression and the surface is the tangent surface to the curve. Therefore, all developable 
surfaces can be represented as planes, cylinders, cones and tangent surfaces, or combinations 
thereof [41]. 

This description of developable surfaces leads to the following duality between space 
curves and developable surfaces [41, p.72]: 


Curves Developable Surfaces 
2 points determine a line. 2 planes determine a line. 
3 points determine a plane. 3 planes determine a point. 
2 consecutive points on a 2 consecutive planes of the 

curve determine a tangent line. family of planes determine 

a characteristic line. 

The curve is the envelope The developable surface is 

of tangents. generated by characteristic lines. 


Gaussian Curvature 
Using equation (4.1) 
R(u,v) = a(u) + vB(u), 


the first and second derivatives of a surface can be represented as 


R, (u,v) = &(u) + vB(u) (4.12) 
R,(u,v) = Blu) (4.13) 
Ruu(u,v) = &(u) + A(u) (4.14) 
Ryy(u,v) = Blu) and (4.15) 
Rika) a0: (4.16) 


Also recall from (2.9) that 
EGS Renan | > 0 


for a regular surface. The surface normal is given by (2.7) 


R, x R, 
Ru x Ry| 
_ (&(u) + vB(u)) x Blu) 
VEG — F? 
_ &(u) x B(u) + o(B(u) x B(u)) 
VEG —- F? 
wes (4.17) 
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N = 


Recall from (2.23) that Gaussian curvature is represented as 


LN — M? 


“= BG= FF 


where, from equations (2.13) and (2.12) 


N = Ryw-N and 
Vi — ete IN 


Thus, from (4.16), MN = 0 and the Gaussian curvature can be rewritten as 


M2 
EG — F2 
= (Roa i N)? 
EG — F2 
(B(u) - &(u) x B(u))? 
(EG — F?)2 
: A)? 
Set 2 ay (4.18) 


K 


Thus, from condition (4.10), zero Gaussian curvature is a sufficient and necessary condition 
for a surface to be developable. 

Since the Gaussian curvature of a developable surface is zero everywhere, the maximum 
and minimum principal curvatures (2.21) and (2.22) of a developable surface can be written 
as 


Kar = 2 + |H], ia el ae 


These principal curvatures reduce to 


Krioe Cl Kee = 0 “when 0. (4.19) 
Rigs = Opwknn = OL | when e=0 and (4.20) 
Kina Ue Kea = coe wien wie) (4.21) 


It is clear from equations (4.19) to (4.21) that at least one of the principal curvatures is 
always zero on a developable surface, which agrees with the fact that the Gaussian curvature 
is zero everywhere (see equation (4.18)). Kmaz and Kmin from (4.19) and (4.21) respectively 
are each termed the nonzero principal curvature, K*, where 


toe 2H. (22) 


Mapping 


Isometric mapping is defined as a mapping that preserves arc lengths between the mapped 
surfaces. ‘T'wo surfaces that can be isometrically mapped to one another are called isomet- 
rics. Since arc lengths are preserved, the first fundamental forms of isometrics are the same. 
Since the Gaussian curvature of a surface and the geodesic curvature of a curve on a sur- 
face depend only upon the coefficients of the first fundamental form, corresponding points 
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on two isometric surfaces have the same Gaussian curvature and corresponding curves on 
isometric surfaces have the same geodesic curvature at corresponding points [26]. 

A conformal mapping is one in which the angle between every pair of intersecting arcs 
is the same on both surfaces. This occurs if the coefficients of the first fundamental form on 
the original surface are proportional to those on the inverse surface [26]. Therefore, every 
isometric surface is conformal. 

Since a plane has zero Gaussian curvature, only a surface with zero Gaussian curvature 
can be mapped to it. Therefore, only a developable surface can be isometrically mapped to 
a plane [26, p.189]. 

Since the tangent planes along any ruling are constant, developable surfaces can be un- 
rolled, or developed, isometrically into a plane without stretching or tearing by successively 
placing each ruling on the plane. One can think of a developable surface as a paper with 
folds. As the number of folds approaches infinity, the surface becomes smooth. The surface 
can then be unfolded again into a plane. This is shown in Figure 4-5. 


Figure 4-5: A Surface with Folds and a Developable Surface 


A geodesic represents the curve with the shortest length between two points on a sur- 
face. Therefore, geodesics on a plane are straight lines. Since developable surfaces map 
isometrically to planes, geodesics on developable surfaces map to straight lines on a plane 
and straight lines on a plane map to geodesics on a developable surface. 


Developable Surface Properties 


The following statements are equivalent necessary and sufficient conditions for a surface to 
be developable. Each has been described above. Thorough proofs can be found in Struik [41] 
and Kreyszig [26]. 


1. Developable surfaces possess the same tangent plane at all points of the same gener- 
ator. 


2. The normal vectors on a developable surface along a ruling are parallel. 
3. A developable surface is the envelope of a one-parameter family of planes. 


4. Gaussian curvature is zero. 
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5. The mapping of a developable surface onto a plane is isometric. 


6. Geodesics on a developable surface map to straight lines on a plane, and straight lines 
on a plane map into geodesics on a developable surface. 
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Chapter 5 


Design With B-Spline Developable 
Surfaces 


5.1 Introduction and Literature Review 


In the current literature, there are three main approaches to representing developable sur- 
faces in Bézier form. These approaches are described below. This chapter explores user- 
friendly methods for the design of developable surfaces. Section 5.2 describes the design of 
developable B-spline strip surfaces that lie between parallel planes. Section 5.3 describes 
the design of developable B-spline surfaces with general three-dimensional directrices. The 
development of a developable surface onto a plane is described in Section 5.4. This is 
an efficient implementation of the method in the current literature, also described below. 
Section 5.5 provides some engineering examples. 


5.1.1 Developable Bézier Strips 


Aumann [1] provides the basis for a user friendly and computationally inexpensive method 
of designing developable surfaces. He represents developable surfaces in terms of two Bézier 
curves (directrices) and rulings between pairs of points from each curve, as discussed in 
Section 4.2, but restricts the two directrices (r4(u),rg(u)) to lie on parallel planes so that 
the tangent vectors r4(u) and rg(u) are parallel, i.e. rg(u) = p(u)ra(u), where p(u) is a 
linear function of u. This condition automatically satisfies the condition (4.9). 


His design philosophy is to consider r,4(u) as a design curve, which means its degree, knot 
vector and control points are specified. In addition, the two end points rg(0) and rg(1) of 
the resulting directrix are provided as boundary conditions. The knots and internal control 
points of the second directrix rg(u) are then computed using a simple system of linear 
equations. He investigates the case when the design curve is a cubic Bézier curve and p(u) 
is a linear function of u. Hence, the resulting directrix is quartic. Figure 5-1 illustrates the 
terminology. 

Frey and Bindschadler {14] extend the work of Aumann [1] by generalizing the degree of 
the directrices and by considering the case where p(u) is quadratic. Aumann [2] further ex- 
tends his work so ra(u) and rg(u) have the same degree. All three papers discuss regularity 
conditions of the developable surfaces. 
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Figure 5-1: A Developable Surface 


5.1.2 Duality Between Points and Planes 


The second approach, introduced by Bodduluri and Ravani [3], is based on the represen- 
tation of developable surfaces as an envelope of a one-parameter set of planes as described 
in Section 4.2. They express the surfaces in terms of plane geometry using the concept of 
duality between points and planes in 3D projective space. They interpret the developable 
surface as the set of its tangent planes represented in the dual form. Their method is math- 
ematically elegant; however, it is somewhat less user friendly and requires a conversion to 
a standard tensor product form [37]. 

Pottmann and Farin [37] generalize the results of Bodduluri and Ravani [3] by providing 
algorithms for converting the dual representation of developables to standard tensor product 
form. Hoschek and Pottmann [23] derive algorithms for interpolation or approximation 
with developable B-spline surfaces. Given a developable surface patch 2, they pick out 
an appropriate number of tangent planes Y; of a developable surface patch 2, and either 
interpolate or approximate this data with a developable B-spline surface. 


5.1.3 Nonlinear Representation 


The third approach is due to Lang and Roschel [29]. They obtain conditions for the control 
nets and weights of rational Bézier surfaces to be developable, which leads to a complicated 
nonlinear system of conditions. Although their work provides a useful criterion for checking 
the developability of a given ruled surface patch, it is difficult to design a developable surface 
with this method. 


5.1.4 Development 


Faux and Pratt [11] provided a method for the development of a curve on a developable 
surface onto a plane based on the fact that a developable surface maps isometrically onto a 
plane. This method is also described by Gurunathan and Dhande [16]. Another method is 
described by Clements and Leon [7], in which a geodesic on the surface is mapped to the 
plane and points on the plane are determined in reference to this geodesic. This method, 
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however, requires the mapping of the geodesic and then the calculation of the location of 
each edge, and is therefore slower than the other method [11]. 


5.2  B-Spline Developable Surfaces with Directrices in Paral- 
lel Planes 


This section further extends Aumann’s work [1] such that n Bézier developable surfaces, 
whose two directrices lie in parallel planes, are strung together by joining them along their 
end rulings. In this manner, they can be represented as a single B-spline surface with C? 
continuity. 

From Section 4.2, a surface is developable if and only if the tangent plane is constant 
along a generator. If the two directrices r4(u) and rg(u) are restricted to lie in parallel 
planes, then the tangent vectors to r4(u) and rg(u) must be parallel at every u for the 
surface to be developable. Therefore, condition (4.9) is equivalent to 


rp(u) = p(u)ra(u) (5.1) 


where p(u) is a scalar function of u. In this thesis, p(w) is selected to be linear, expressed 
as 


p(u) = po(1 — u) + piu. 


The design philosophy is to consider r4(u) as a design curve [1, 14], which means its 
degree, knot vector and control points are specified. Also, the two end points rg(0) and 
rg(1) are provided as boundary conditions. The knots and internal control points of the 
second directrix, rg(u), are then computed. Since multiplication of B-spline functions is 
complex and computationally expensive, the Bézier segments are first extracted from the 
B-spline curve by knot insertion using Boehm’s algorithm [22] described in Section 3.3. 

Figure 5-2(a) shows the design curve r4(u) which is a planar cubic B-spline curve with 
knot vector T = (0000 4 2111 1) and control points (0, 2, 0), (0.8, 0.5, 0), (1.6, 1.5, 
QO), (2.4, 1.5, 0), (2.6, -0.5, 0) and (4, 0, 0). Figure 5-2(b) shows three Bézier segments 
after knot insertion with the two endpoints of rg(u), which are supplied by the user to be 
r3(0) = (—1,3, —3) and rg(1) = (2, 1.3, —3). Note that r4(u) is in the z = 0 plane and the 
rg endpoints are in the z = 3 plane, parallel to the plane of r4. 

Let r4(u) be a degree n curve with m Bézier segments. Since p(w) is linear, the degree 
of rp(u) isn +1. We denote the jth Bézier segments as 


n 
r,(u) = S- ajn+iBin(u) 7=Uee.,m—1 and 
i—0 
n+l 
rau) = d dinstiBinti(u) j=0,...,.m—1. (5.2) 
— 


The control points a and b are two-dimensional vectors since they have been restricted to 
lie in parallel z planes. Substituting (5.2) into equation (5.1) yields 


nr 


(n+1)) (Bineysit1 — byn4i4:)Bin(u) 
7=0 
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(a) Design Curve r, (u) with Contro! Polygon and Knots (b) Design Curve after Knot Insertion with Bezier Segments 
and Endpoints of r, (u) Shown 





(c) Design Curve and Resulting Directnix before Knot Removal (d) Final Developable B-Spline Strip Surface 


Figure 5-2: Design of Developable B-Spline Strip Surface 
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7=0 
n—-1l 
= n> (po(l—u) + pit) (ajnsiti — Ajngi)Bin-1(u) (5.3) 
2=0 


for 7 = 0,...,m—1. Using the definition of a Bernstein polynomial (3.1), 
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Bi+in{u). (5.5) 


Substituting (5.4) and (5.5) into (5.3) yields 


(n ts 1) S (BD j¢n41)4i41 a b5(n41)+i) Bin (u) 
1=0 


n—1 
= n 3 [po(1 = 1b) eet) ae p\uBin—-1(u)] (Ajn+i+1 7 Ajn+i) 
1=—0 


n—l 
= yy [po(n — 7) Bin(u) + pilt + 1) Bisin(u)] (ajntiti — ajnti) 
i= 
= por(ajn+1 os ajn)Bon(u) 
n—1 
ap ye [Po(n — t)(@jnsit1 — Ajnti) + pit(Ajnsi — ae Bin(u) 
2 
+pin(ajnin at Ajn+n—1)Pnn(u)- 
Since Bernstein basis polynomials are linearly independent, each set of coefficients must be 
equal, yielding the following system of n — 1 equations [14] for the jth Bézier segment 
npy(Ajnti —Ajn) = (n+ 1)(Din4141—Pyn1)) (5-6) 
(n — i) po (Ajn+i+1 — ajn+i) + 1p} (Ajnti —Aajnzi-1) = (n+ 1) (Dj(n41)4i41 7 bj(n+1)+i) 
np (again — AGg4iyn—-1) = (2 +1) (BG 41y(ngay — bG4ry~mgy—1)s 
where 2 =1,...,n—1. 


Now consider joining two consecutive Bézier patches with C? continuity. First, C! 
continuity requires that the first derivatives along the connecting edges of both patches 
must be equal. Since the connecting edge is a ruling and is therefore linear, this requirement 
reduces to 


i eer) (5.7) 
(1) = 210). (5.8) 
From equation (5.1) we have 
PB(1) = o(1)FA(1) 
te (0) = po Cyrt*(0). (5.9) 


Substituting equations (5.9) into equation (5.8) and using equation (5.7) yields 
po) = pitt. (5.10) 


ol 


C? continuity requires that the second derivatives along the connecting edges must be 
equal. Again, the connecting edge is linear, which reduces this requirement to 


Py = Ra) (5.11) 
Pe(1) = #R°'(0). (5.12) 
Taking the derivative of equation (5.1) yields 
Pp (u) = p(u)i4(u) + o(u)éa(u) (5.13) 
where 
p(u) = pi — po. 


Evaluating equation (5.13) at the connecting edge yields 


(1) = pit (1) + (ol - 08)2,(1) 
) 


Substituting equations (5.14) into equation (5.12) and using equation (5.11) yields 
= tp (5.15) 


Each Bézier patch has n+ 1 control points along r4(u) and n + 2 control points along 
rg(u), where n is the degree of the original B-spline design curve. The system of equations 
for a single patch derived in equation (5.6) can be solved for the coefficients of the second 
directrix 





Dj(n41)41 _ Dj(n+1) + ae PA Ajn+1 — ajn) 
—1 
Dj(n41)4i41 Dy(n4iyti + aa Po (Ajnsit 7 Ajn+i) 
a 
+—— pi (@jnti — Ajnti- ar t=1,...,n-1 
by+i(nti) = Bytiynen-1 + — Play tayn — Aj 41)n-1) (5.16) 


where 7 = 0,...,m—1 for a system of m patches. Adding equations (5.16) for m patches 
yields a vector equation that can be solved for p§ and Dees in terms of the first and last 
points of the second directrix: 


Dm(n+1) —bo = Povo ar maa (32%) 


where 


vo = wey est Eee) and 
ton - waar (3+ 2) 


o2 


V1 


Combining equations (5.10) and (5.15), we can express p and pitt in terms of the two 
end coefficients p§ and p7"~' as follows: 


pt=g=— po 7 =0,...,.m—1. (5.18) 


J AO 4 





Therefore, given a design directrix, r4(u), and the endpoints of the second directrix, rg(0) 
and rpz(1), all the interior control points of the second directrix can be obtained succes- 
sively using equations (5.16). After the second directrix is calculated, the first directrix is 
degree elevated so the two curves may be combined into a B-spline strip. Any knots that 
may be removed from the surface while leaving the curve geometrically and parametrically 
unchanged are removed using the knot removal technique [42] described in Section 3.2. 

Figure 5-2(c) shows the computed rg(u) together with r4(u) in Bézier form and the two 
end rulings. Figure 5-2(d) shows the final B-spline developable surface after knot removal. 

Developable surfaces become irregular when |R, x R,| vanishes, which leads to the 
regularity condition for a strip surface [14] 


p(u) > 0, Oi eee 


From a user’s point of view, once he/she determines one of the endpoints, it is convenient 
to know the range that the other end point can take. From (5.18) it can be seen that 
if both p§ and an are positive, then all p coefficients are positive and the surface is 
regular. The coon p can be determined by solving the vector equation (5.17). The 
resulting equations for p? and p7’~ ' must be positive for the surface to be regular, yielding 
the following restrictions: 

(OF (n+1) — bj)vi - (Br in+1) — b9)ot ee 


0 
; Upri — YIU 


Po = (5.19) 


and 


(0 ity — 06)V6 — (OZ ns) — 58) UG 
tna mine) ~ Po) = Ominsa) — *a)¥0 =) (5.20) 


VOY] — VT UO 
where the superscripts denote the x and y components of the vectors. Given one endpoint, 
a regular surface may be formed by selecting the second endpoint to lie within the range of 
values determined by equations (5.19) and (5.20). 

It is also interesting to observe [14] that when the two vectors joining two end points of 
each directrix are parallel and their magnitudes are multiples of a constant, the resulting 
developable surface is a generalized cone. When the constant is equal to one, the surface is 
a cylindrical surface. In other words, if p(u) is a constant not equal to one, the surface is a 
generalized cone, and if p(u) = 1, the surface is a cylinder. 


Algorithm 


The following steps describe the algorithm for the solution method. 
1. The user provides the design curve r4(u) and the endpoint rg(0). The program 


provides a range of values for rg(1) that ensure a regular surface, and the user supplies 
a specific value for rg(1). 


o3 


2. Bézier segments are extracted from r4(u) by knot insertion. 
3. The system of equations (5.16) is solved for the internal control points of rg(u). 


4. Once the resultant directrix rg(u) is computed, the developable Bézier surface is 
obtained by degree raising r4(u) and substituting in equation (4.2). 


9. Knots that can be removed while leaving the surface unchanged are removed using 
the method described in Section 3.2. 


5.3 Developable Surfaces with 3D Directrices 


5.3.1 Developable (n,1) Bézier Surfaces 


In Section 5.2, we restricted the two directrices to lie on parallel planes. In such case, r4(u) 
and rg(u) must be parallel to satisfy condition (4.9). In this section, the directrices are not 
constrained to lie on parallel planes, so equation (5.1), rg = p(u)ra, no longer applies. 

The design philosophy here is the same as that of the previous section. Namely, r4(u) 
is designated the design curve with its degree, knot vector and control points specified. The 
two end points rg(0) and rg(1) are provided as boundary conditions and the knots and 
internal control points of the second directrix, rg(u) are computed. 

We begin this development by assuming the design curve r,4(u) is a Bézier curve of 
degree n, 


n 
r4(u) = aera), CSS 
=5 


where a; = (a7, a¥,a?)? are control points. The resultant curve rg(u) is also a Bézier curve 


of degree n. In Section 5.3.2 we investigate the case when the curves r4(u) and rg(u) are 
B-spline curves. If b; = Ce 1 = 0,...,n, are the control points of the degree n 
Bézier curve rg(u), and a;, b; and c; are 


= hae ie ; 

a; = "dene —a7—a, i=0,...j72—1 

im LD tie . 

b= (b7, bY, bF) —Ipeee— bp; 1— 0,...,%— 1 
a : 

Co = (ee aes — Ib, a; dale 


then equation (4.9) can be expressed as 


n n—-1 n—-1 
(>, e.Bin() x (1 SS 8.Bn-a(u) ( yo BiB, n-1(u) =al); (522) 
1=0 I=) 7=0 


Since equation (5.21) is a scalar polynomial equation of degree 3n — 2, we can rewrite it 
(upon division by n”) in the form 


3n—2 
>, Ex Bran-2(u) = 0 (5.22) 
k=0 


o4 


where 


min(2n—1,k) Ga (ea 


Ey, = » (P?bz_; - Pb, 7 P?bi_;) (5.23) 
i=maz(0,k—n+1) ( k ) 
and 

“a Ov) 

Fit > (4%-5- 395-3)“ an=ty 
jg=marz(0,i—n+1) 1 

io see LZ Z—2 (3) (773) 

Py = x (C5 G;_5 — 591-3) 
j=marz(0,1—n+1) 1 

min(n,t) @ Ea 

—— 3 (cra!_, ~ at_,;)-4— 


jg=maz(0,1—n+1) 


using the algorithms for multiplication of polynomials in Bernstein form described by 
Farouki and Rajan [10]. Since equation (5.22) is valid for all u and the Bernstein basis 
functions are linearly independent, £; must be zero for all k, 0 < k < 8n—2. This leads us 
to a system of 3n—1 polynomial equations with 3(n—1) unknowns, i.e. b; (¢ = 1,...,n—1) 
where the two boundary points, bo and b,, are given. Therefore, there are two more equa- 
tions than unknowns and the system is overdetermined. The system of equations was solved 
using a Numerical Algorithms Group program [35] for minimization. 


Solution Method Using Unconstrained Optimization 


One way to solve the overdetermined system is to seek a least squares solution to E(x) = 0 
where 


bY 


n—1? 


eee (DOs 0 De 


z 
n—1l? b 


n—l 


ie 
and 


E(x) = (E(x), Ey (x), oe , A3n-2 (x))?. 


This can be obtained by minimizing the function 


F(x) = = (ER (x) + B3(x) +... + £3, o(x)) = 5B" (x)E 


orn 


aly , (5.24) 
7=() 


The problem can then be reduced to searching for zeros of the gradient vector field V F(x), 
1.e. 


VF (x) = 0. (5.25) 


The Newton-Raphson method can be used to iteratively solve the nonlinear system (5.25). 
The Taylor series expansion about the current point is 


V F(x + 6x) = g(x) + [H(x)]6x (5.26) 


DO 


where g(x) and {H(x)] are the gradient vector and the Hessian matrix of the objective 
function F(x). Using (5.26), the step direction 6x is found from 


[H(x)]éx = —g(x): (5:27) 
The value for the subsequent iteration, x is 
x=x+ 6x. 


Using the special properties of a least squares function such as equation (5.24), the deriva- 
tives of F(x) are represented as 


g(x) = VF (x) = [A(x)]E(x) 


and 


3n—2 
[H(x)] = V7F(x) = [A(x)][A(x)]? + $°> E,(x)V7E;(x) 
0 


where 
CA) = ey eee | 


is the (3n — 3) x (3n — 1) Jacobian matrix, the columns of which are the first derivative 
vectors VE; of the components of E ({A;;] = 52). 
In any minimization process, the initial approximation for x greatly affects the resulting 


solution. Here we use a linear interpolation to obtain the initial values for b;, where 


Op a (Ope an eee (oe) ten rit a 


Simple Bounds 


In many applications, it is helpful to introduce simple bounds on the variables [15] such 
that 1; < x; < u;, where J; and u; are lower and upper bounds, respectively. This is a 
constrained optimization problem; however, it has a particularly simple form. A solution 
x* is termed feasible if it satisfies all constraints. A step 6x is in a feasible direction if it 
does not cause x to exceed any of the bounds. 

The following discussion treats a generic iteration, so the iteration number has been left 
out for simplicity. At any iteration, some of the variables will be equal to their upper or 
lower bounds. These variables xpx, are fixed on their bounds and are consequently termed 
fized variables. The associated bounds are placed in the active set. The remaining free 
variables Xp are allowed to vary within the bounds. The search direction 6x is chosen to 
change only the free variables. Therefore, 6xpx = 0 and equation (5.27) reduces to 


[He_a(x)|éxer = —rr(x). (5.28) 


With dxpp we update x for the subsequent iteration as 


auf am hie 


Xrx 
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If dzeR forces Zpp to violate one of the bounds, that bound is added to the active set for 
the subsequent iteration. 

At each iteration, the active set must be analyzed to determine if the constraints in the 
active set should remain active. A solution x* is a local minimum only if there is no feasible 
step dx that results in a lower value. Thus, if 


(6x)? g(x) (5.29) 
is non-zero for any 6x, then x* is not a local minimum. Then (5.29) can be represented as 
g(x).= [A]’A 


where [A] is the set of constraint coefficients and \ are Lagrange multipliers. For simple 
bounds, [A] is composed of 1 for J; and —1 for u;, so the Lagrange multipliers can be 
estimated as 


_ Prix for 7: — 1, 
ae —gi(x) forz;=u; — 30) 


Thus, if a Lagrange multiplier for an active simple bound is negative, the corresponding 
bound should be removed from the active set for the subsequent iteration. 





(a) Unconstrained Cylindrical Developable Surtace (b) Cylindrical Developable Surface 
With Bounds 


Figure 5-3: Effect of Simple Bounds on Developable Surfaces 


Many solutions are possible for a single set of equations. For example, Figure 5-3 shows 
two equally viable solutions to the same set of equations. Both examples have design curves 
with control points (0, 0, 0), (1.8, 3, -0.7), (3.3, -2, 1.5) and (4, 0.5, 0) and endpoints (0, 0, 
3) and (4, 0.5, 3) for the second directrix. Figure 5-3(a) is unconstrained, yielding interior 
control points for rg of (1.8, 3, 3) and (3.3, -2, 3). Figure 5-3(b) is unconstrained in the z 
and y directions, but has a lower bound equal to 5 set in the z direction. This produces a 
resultant curve with interior control points (1.8, 3, 5) and (3.3, -2, 5). Since both cylinders 
are developable surfaces with C? continuity that include the endpoints specified by the user, 
they are equally valid solutions to the equations. Both surfaces are exactly developable with 
zero Gaussian curvature. Although this example is rather trivial, it demonstrates the control 
available to the user to affect the shape of the resulting surface by specifying bounds on the 


of 


solution. Another example is given in Section 5.3.2. 


Generalized Elimination Method 


From equation (5.23) it can be seen that Eo(x) = 0 reduces to a linear equation, 
(bo — ap) x (a, — ap) - (bi — bo) = 9, (5.31) 


where b, is the only unknown. This condition forces the first four control points, namely 

aj, a1, bo and bj, to lie on a plane, and has the following geometric interpretation. The 

two Bézier curves r4(u) and rg(u) pass through the points ag and bg and their tangents 

are in the directions of the vectors aga, and bob,. Since a developable surface has the same 

tangent plane at all points along a ruling and these tangents are on the same ruling, u = 0, 

they must be coplanar. Therefore, aj, a;, bg and b; must lie on the same plane. 
Similarly, the last equation, £3,-2(x) = 0, is linear in by-): 


(b, — ay) X (apeeanet) 40, — bni1) =0. (5.32) 


The linear equations (5.31) and (5.32) can be used to reduce the dimensionality of the 
optimization by two. The new objective function can be introduced without terms E2 (x) 
angels a5 oc)a 1.6. 


F(x) = E?(x) + Ef(x) +... + E%,_3(x). 


The following discusses the generalized elimination method [13, 15], in which the number of 
variables in the optimization is reduced by the number of independent linear constraints. 
In general, if we denote the number of variables in the optimization as N (ie. x € RY) 
and the number of linear equations as M, the set of linear equations can be written as 


ae 


in which [A] is M x N and b € R™. In other words, the ith row of [A] contains the 
coefficients corresponding to the zth linear equation, and is therefore the gradient of that 
equation. 

This reduction in dimensionality can be described formally in terms of two subspaces, 
which can be described as 


e the M-dimensional subspace defined by the rows of [A]. 


e the N — M complementary subspace of vectors orthogonal to the rows of [A]. 

Let [Y] denote any N x M matrix whose columns form a basis for the range space of [A]. 
In fact, [Y] may be taken as [A]?. Let [Z] denote an N x (N — M) matrix whose columns 
form a basis for the set of vectors orthogonal to the rows of [A], so that [A][Z] = [0]. 
The matrix [Z] has linearly independent columns 2), Z2,..., Z2N—m which are in the null 
space and act as basis vectors for the null space. Since [Y] and [Z] define complementary 
subspaces, every N-vector x has a unique expansion as a linear combination of the columns 
of [Y] and [Z] 


x = [Y]x, + [Z]x., (5:33) 
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where M-vector x, is the range-space portion of x and (N — M)-vector x, is the null-space 
portion of x. By using the chain rule on the derivative of equation (5.33), it follows that 


Ve, = (Z]’ Ve. 
Thus, 
Vz, F(x + [Z]x) = [Z]’ V2F (x) = [Z]’ g(x) 
is the reduced gradient vector and 
Vz, F(x + [Z]x) = [Z]? VF (x)[Z] = [Z]’ (A(x) [Z] 


is the reduced Hessian matrix. Given a feasible initial point x9 we solve a reduced size 
(N — M) x (N — M) linear system 


(Z]" [H(x)|[Z]x. = -[Z]’ g(x) 
for x,, where the subsequent iteration is 


mee Albee 


5.3.2 B-Spline Developable Surfaces 


In this section we merge m Bézier patches together by joining them along their end rulings 
with C? continuity such that the resulting developable surface can be represented by a single 
B-spline surface. The design procedure is the same as that of the Bézier case, except that 
the design curve is a B-spline curve. Since multiplication of B-spline functions is complex 
and computationally expensive, we first extract the Bézier segments from the B-spline curve 
by knot insertion [5] as described in Section 3.3. We then evaluate equation (5.23) for each 
Bézier segment and add the resulting equations to the objective function after squaring 
them. If equations Ep and £,,(3n-2) are used in the general elimination method, they are 
not included in the objective function. 


Recall from Section 3.2 that, for a degree (n, 1) B-spline surface with m Bézier segments, 
the C! and C? continuity conditions at the joints are given by 


hilee ei ba:) — henry — Wei), eae? = ne (5.34) 
and 
Ais h; 
bri-1 + ——(Dni-1 — Pai-2) = Pnigi + — (Bnidi — Dri+2). (5.35) 
ie Qo 


Equations (5.34) and (5.35) reduce to six scalar linear equations for each junction. Using 
these linear equations in the generalized elimination method described in Section 5.3.1 
ensures that the surfaces will be exactly C? continuous; otherwise they may have some 
residual value after optimization that would cause discontinuity. For a degree (n, 1) B-spline 
patch with m Bézier segments, the number of unknowns and linear equations is 3(mn — 1) 
and 6(m — 1) + 2 respectively, so the reduced number of unknowns is 3m(n — 2) + 1. 


For a B-spline surface, the initial approximation for x is calculated in the same manner 
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as for the Bézier case, expressed as 


bj = a; + ——(bo — a9) + ———(bmn — Amn), §=1,...,.mn—1. (5.36) 
mn mn 


Once the control points of the second directrix are computed, any knots that may be 
removed from the curve while leaving the curve geometrically and parametrically unchanged 
are removed using the knot removal technique described in Section 3.3. 

Figure 5-4(a) shows the design curve r4(u) which is a cubic B-spline curve with control 
‘points (0, 4, 0), (1, 2, 1), (2, 4, 1), (3, 1, 0) and (4, 2, 0) and knot vector (0, 0, 0, 0, .5, 1, 
1, 1, 1). Figure 5-4(b) shows the design curve after knot insertion with the resulting two 
Bézier segments and the selected endpoints rg(0) and rg(1), which are (0, 4, 3) and (8, 0, 3). 
When this problem is treated as an unconstrained optimization, the answer converges to the 
trivial solution where the internal control points of the new curve are equal to the internal 
control points of the design curve as shown in Figure 5-4(c). To avoid this situation, lower 
bounds are set on the z variables where /; = 3.0, 7 = 2, 5, 8, 11 and 14, and the resulting 
internal control points are (2, 0, 5), (4, 4, 5) and (6, -2, 3). The resulting surface is shown 
in Figure 5-4(d) with two Bézier patches and the control mesh delineated. Figure 5-4(e) 
shows the final developable B-spline surface after knot removal. 

Aumann [2] proves that a degree (3-1) Bézier surface whose directrices lie in parallel 
planes is developable if and only if one of the following conditions holds: 


1. The rulings R(0,v) and R(1,v) are intersecting or parallel. 
2. The coefficient of the cubic term of the design curve is zero. 
3. The design curve contains exactly one singular point. 


In constructing examples for this thesis, it was found that skew end rulings induce non-zero 
Gaussian curvature into (3-1) Bézier surfaces whose directrices are 3D space curves as well. 


Algorithm 


The following steps describe the algorithm for the solution method. 


1. The design curve.r4(u) and the endpoints rg(0) and rg(1) are provided by the user 
along with any constraints such as simple bounds or constraining planes on the interior 
control points of rg(u). 


2. Bézier segments are subdivided out of the B-spline curve by knot insertion as described 
in Section 3.3. 


3. Initial values for x = (bf, 66, 0§,---, btn, Yin, O71) are calculated using a linear inter- 
polation method (5.36). 


4. The objective function is optimized in two steps. First, beginning with the initial 
values, an initial feasible point is calculated that meets all constraints and bounds. 
Then, an optimal solution is found following the method described in Section 5.3.1, 
using the following steps: 


e If x satisfies the objective function F(x) = O and the Lagrangian multiplier 
estimates are non-negative for all active bounds, the iteration process is complete. 
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{a) Design Curve ¢ a(u) with Control Potygon and Knots (b) Design Curve after Knot Insertion with Bezier Segments 
and Endpoints of r g(u) Shown 





(c) Resulting Surface Without Simple Bounds 





a " . a a 


(d) Design Curve and Resulting Directrix before Knot Removal (8) Final Developable (3,1) B-Spiine Surface 


Figure 5-4: Design of Developable B-Spline Surfaces 
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(Note that the objective function is a sum of squares, so the minimum possible 
sum is zero.) 

e A feasible step direction 6x is found using (5.28). 

e If 6x causes x to exceed any of the bounds, those bounds are added to the active 
set. 

e The Lagrangian multiplier estimates are calculated for all bounds in the active 
set using (5.30). If any are negative, the associated bound is removed from the 
active set. 


e x is set to x + 6x and the next iteration is begun. 


5. Once the resultant directrix rg(u) is computed, the developable B-spline surface is 
found using equation (4.2). 


6. Knots that can be removed while leaving the surface unchanged are removed using 
the method described in Section 3.2. 


7. The Gaussian curvature bound is calculated using the method described in Sec- 
tion 5.3.3 to gauge the quality of the resulting surface. 


8. The surface is developed onto a plane using the technique described in Section 5.4. 


5.3.3 Accuracy of Resulting Developable Surfaces 


Since the optimization results do not always precisely satisfy F(x) = 0, we need a measure 
to determine how close the resulting surface is to an exactly developable surface. Recall 
from Section 4.2 that a necessary and sufficient condition for a surface to be developable is 
for the Gaussian curvature to be zero. This property can be used as a measure. For each 
resulting Bézier patch R(u, v) = S-i~o Yj=0 R,; Bin(u)B;,1(v), the Gaussian curvature, K, 
can be expressed in terms of ratios of bivariate Bernstein polynomials as [31] 


=(Se : 
K (uy) = S-Bas)(S Rev) ~ (8 “Rae 


Gre v) 
Keno v) 
ae eee Ky; Bi en—4(u) Bya(v) 


= (5.37) 
Sy DD ro een aen) 3 aa) 


where S = R, x R, and S = |R, x R,|. If R(u,v) is a developable surface, then K;; must 
be zero fori = 0,...,6n —4 and 7 = 0,...,2. However, for surfaces where F(x) 4 0, the 
Kj; coefficients are not necessarily always precisely zero. 

From equation (5.37), we have 


rand) 
K denom (u, v) 


— |Knum(u, v)| 


— \Kdenom(u, v)| . Ce 


JK (u,v)| = 





We can easily find the upper bound of |Knum(u, v)| as 


6n—-4 2 


ee 6G bron) 


7=0 7=0 


|Knum(u,v)| = 
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6n—4 2 


< >” YS  |Kij| Bien-4(u)Bj,2(v) 
2=0 7=0 
6n—4 2 
maa. 3 S— Bisn—4(u)Bj,2(v) 
v 7=0 7=0 
= max |Kj)|. 
tJ 


The denominator Kgenom(u,v) is always positive if the surface is regular, since |R, x R,| 
is always positive for a regular surface. However, some of the coefficients 5;; can be zero 
or negative. If there are non-positive S;; we can subdivide the surface into smaller patches 
until all S;; become positive. This is always possible by virtue of the convergence of the 
control polyhedron to the surface with repeated subdivision [22]. Then if all S;; > 0 we 
have the inequality | 


8n—-4 4 
|Kaenom(, v)| — Se S— Siz Bign—4(u) B;,4(v) 
t=O 70 
8n-4 4 
> min|S;;| S- S_ Bign—4(u)B;,4(v) 
2 iO —0 
— Seino 
ee 
Therefore we can rewrite equation (5.38) as 
Pen 
Tay ea aa (5.39) 
min;; 5;; 


If min;; S;; is very small, the bound of the Gaussian curvature can be extremely large. 
In such case we further subdivide the patch and re-evaluate the bound until such bound 
converges to a stable value. 


Alternatively, the curvature may be analyzed at discrete locations on the surface, which 
provides a rough estimate but may not find the actual bound of the maximum curvature. 


We use the Gaussian curvature bound to give a physical measure of the developability of 
the surface. Throughout this thesis we say the surface is exactly developable if |K| < 10719. 
For example, the curvature of the surface shown in Figure 5-4, as determined using ten sub- 
patches per Bézier patch in equation (5.39), is zero, so the surface is exactly developable. 


5.3.4 Developable Surfaces Constrained by Planes 


In many engineering applications, one or both of the directrices are restricted to planar 
curves which are not necessarily parallel. In this situation, we can force the control points 
b; and thus the directrix rg(u) to lie in a plane. This plane is determined by the two end 
points, bp and bmn for a degree (n,1) B-spline patch with m Bézier segments, and a point 
g chosen by the user. The equation of the plane is given by 


az+by+cz+d=0 (5.40) 
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where, for a degree (n, 1) B-spline patch with m Bézier segments, 


(6 Si gq? (be a g°) cae (b6 - Coon = g°) 
(69 — 9°) (Bin — 9°) — (65 — 9° )Omn — 9°) 
(b6 _ g Hoon a g’) ~ (b5 ~ GOs _ g°) 
—ag’ — bg’ — cg’. 


QS Gears 
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For each control point b; to lie on the plane, its z, y and z components must satisfy 
equation (5.40). These linear constraints can be used to reduce the number of unknowns in 
the equation using the generalized elimination method described in Section 5.3.1. 

However, if the (mn—1) planar constraints, the (6m-—6) C! and C? continuity conditions 
and the two linear developability conditions Eo and £,,(3n-2) are all used to reduce the 
number of variables, the problem may become so overconstrained that is is difficult or 
impossible to find a feasible solution. For example, a cubic B-spline design curve with 
two Bézier segments has 15 unknowns (bj), ba, b3, b4, bs) and 13 linear equations, leaving 
just 2 variables for the entire problem. In this case, we use a subset of the total available 
linear equations as linear constraints and add the remaining equations to the objective 
function after squaring them. In selecting the linear equations to be used as constraints, 
it is important to remember that the linear constraints may be satisfied exactly while 
the objective function equations may not be. For example, if C? continuity is a must, 
equations (5.34) and (5.35) are used as linear equality constraints and the equations for 
planar constraints are included in the objective function. 





(a) Developable Surface with Constraining Plane (b) Final Degree (3,1) Developable Surface 


Figure 5-5: Planar Constraint on Developable B-Spline Surface Design 


Figure 5-5 shows a developable B-spline surface that was designed with the second 
directrix constrained to a plane. The same design curve and endpoints were used as those 
used for Figure 5-4, except the second directrix was constrained to lie in the z = 3 plane. 
Note that since the second directrix was constrained to a plane, the trivial solution was not 
an option, and simple bounds were not required. The continuity conditions were used.as 
linear constraints. 


5.3.5 Degenerate Developable Surfaces 


Sometimes topologically triangular patches are required to produce a desired surface design. 
A three sided patch parameterized over a rectangular domain can be constructed by allowing 
one of the boundary generators to be zero length, say R(0,v). Faux and Pratt [11] showed 
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that the surface normal at this point may be evaluated using a limiting process based on a 
Taylor series expansion 


R,(0,v) x Ruy (0,v) 
IR(0,v) x Ruy (0, v)| 


For a degree (n-1) developable Bézier patch, equation (5.41) reduces to 


hen) (5.41) 
u—0 


(P10 — Poo) x (Pi — Pio) 


N(0,v) = 
oe) (ero bone (ii — Pio) 


Therefore, the unit normal vector at the degenerate point always exists as long as the three 
control points, Poo, Pi,0 and P,,1, do not coincide with each other, and it is orthogonal to 
the plane formed by the three control points as shown in Figure 5-6. 


The Gaussian curvature at the degenerate point cannot be calculated using the method 
described in Section 5.3.3 because the denominator, |R, x R, |, vanishes. Wolter and 
Tuohy [43] devised a method for determining Gaussian curvature at a degenerate point using 
the second order partial derivatives of a local height function representation of the surface 
assuming the surface has a tangent plane and well defined curvatures at the degenerate 
point. For this description, assume that the degenerate point occurs at R(0,v). The 
coordinate system is shifted such that the origin coincides with the degenerate point, the 
z axis is aligned with the normal to the surface and the z axis is aligned with the tangent 
to the surface in the non-degenerate parameter direction; in this case, R,(0,0). The y 
axis is orthogonal to the others to form a right-handed system. Three linearly independent 
isoparameter curves in wu are selected on the surface such that a,(u) = R(u,vp), n = 


1,...,3. The normal curvature in the direction 
OR(0,un) _ Oap,(0) 
Wn = — a = ae = (fasua) 


at the point R(0,0) is given as 


n-k, = n-a‘(0) 


Grn 0)x2 3 ogr ANG; Nan =I Coie 


which can be represented in matrix form as 


[L]a=b (5.42) 
where 
a = (Gras 0cusGiy ae 
b = (n-a/(0),n-a5(0),n- a3(0))* 
and 
at Qniy y? 
[L]= | 23 2zey2 y5 


az 22343 3 


If the fundamental forms are converted to the new coordinate system, the Gaussian curva- 
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ture can be represented as 


K = grrQyy — Gen (5.43) 


The Gaussian curvature at a degenerate point can thus be determined by solving (5.42) for 
a and substituting the resulting values for grr, gry and gy, into (5.43). 

To estimate the maximum absolute value of the Gaussian curvature bound of a degen- 
erate surface with the degenerate point at R(0,0), the surface is first subdivided at u = e, 
where € is a small positive number so that the resulting patch with e < u < 1 1s a nonde- 
generate patch. The Gaussian curvature is calculated at u = 0 using equation (5.43), and 
the bound for the Gaussian curvature for the nondegenerate patch is evaluated using the 
technique described in Section 5.3.3. 





Figure 5-6: Triangular Degenerate Cubic Bézier Patch 


Figure 5-6 shows a degenerate cubic-linear Bézier patch with non-planar directrices. The 
design curve is made up of the control points (0, 0, 0), (1.8, 3, -0.7), (3.3, -2, 1.5) and (4, 
0.5, 0). The selected end points are (0, 0, 0) and (4, 0.5, 2). The Gaussian curvature at the 
degenerate point was found to be 8.91 x 10~!*. The surface was split at the parameter value 
u = 0.0001 resulting in a maximum bound for the Gaussian curvature of zero. Therefore, 
the surface is exactly developable. Note that there are only three points which define the 
end rulings and they always lie on the same plane. Therefore, there is no skew between the 
two end rulings and the surface often results in an exact developable surface. 


5.4 Development of a Developable Surface onto a Plane 


In the manufacture of developable surfaces, it is necessary to calculate the plane development 
of these surfaces. This determines cut information for the flat piece of metal so it can be 
cut to the proper dimensions before it is rolled. 

This formulation follows the method of Faux and Pratt [11, p.274], using the fact that a 
developable surface may be mapped isometrically onto a plane. The Frenet-Serret formulae 
state that 


(=A (5.44) 


66 


where t is the unit tangent vector to a curve, n is the unit normal vector to a curve and 
« is the curvature. For a planar curve in the (X,Y) plane, we can define the unit normal 


vector as 


me x t (5.45) 


where e, = (0,0,1). Substituting (5.45) into (5.44) yields 


d2X dY 
i 5.46 
ds? * ds ( ) 
d7Y dX 
—— —_ = 5Ay 
ds? ca ds f ( ) 


where (X, Y) denote the 2D coordinates on the developed plane (X,Y). If we rewrite equa- 
tion (5.46) in terms of the parameter u, we obtain 


Oe de Bd dX od, 
dst ds ao] 
_ d (dXdu\  _dY du 
- a ee 


a 
~ du? \ds du \ds?} “duds 
If (du/ds)* 4 0, this reduces to 
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—- —_- — = 0. 5.48 
du? i can du (4) du (5.48) 
ds ds 
Similarly, equation (5.47) can be represented as 
ry (tx 
iiae (4, ) ay na = (5.49) 


du? * (#) a (#) ee” 

ds ds 
The development is based on the fact that curves on isometric surfaces have the same 
geodesic curvatures. It can be seen from Figure 2-4 that for a planar curve, since the angle 
between the surface normal and the normal to a curve is 90 degrees, the curvature vector 
equals the geodesic curvature vector. Therefore, « in equations (5.48) and (5.49) can be 
replaced by Kg, the geodesic curvature of the curve on the developable surface. If we choose 
the curve on the developable surface to be an isoparametric curve in terms of u, R(w, vp), 
we can replace du/ds and d?u/ds? in equations (5.48) and (5.49) by 


du 
ds 
d*u 
ds? 


1 
Richie 
_ Riu, vn) Ruu(u, Un) 


Ri(u, oe) ; Ri(u, Un) 
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Thus we have 


dX _ 
du ee 
ay _ 
du ae 


dp (Ry . Ruu) 

SES Eo, | 

du ~ (RR) “*9/Pe 

dq (R, ; Ryu) 

— a 5.90 

du qd (R,, Ty) + pkg|Ry| ( ) 

The surface will be developed onto the plane oriented such that the point R(0,0) is 

located at the origin of the (X,Y) coordinate system in the plane and the vector R,,(0, 0) 
coincides with the X axis, as shown in Figure 5-7. We integrate the system (5.50) along 
the directrix that corresponds to v = 0 from u = 0 to u = 1, using, for example, a Runge 
Kutta scheme. Since isometric maps are conformal, the angle between the directrix and the 
generator at (0,0) is the same in both representations and can be found by 


Ru(0,0)  R,(0,0) 


98? =TR.,(0,0) | 1R4(0,0)] 

Straight lines on the surface map to straight lines on the plane, so the generator length 
on the surface is measured and transferred to the plane, giving the initial point for the de- 
velopment of the second directrix. Then the system of ordinary differential equations (5.50) 
is integrated using a fourth order Runge-Kutta method [38] along the directrix that cor- 
responds to v = 1 from u = 0 tou = 1. Finally, the endpoints of the two directrices are 
connected. 

To verify that the development of the surface is accurate, one can check the results 
by computing the angle between the directrix corresponding to v = O and the generator 
corresponding to u = 1, computing the length of that generator, and determining the 
location of the point R(1,1) in the plane. Comparing this point to the point calculated 
above determines the accuracy of the development. The development will be more accurate 
with increased steps in the Runge-Kutta scheme. 

The surface that was designed in Section 5.3.2 is developed onto a plane here. This 
surface is shown in Figure 5-7(a). It was developed onto a plane using 31 steps in a fourth 
order Runge-Kutta routine as shown in Figure 5-7(b). The accuracy for the check at R(1, 1) 
was 107°. 


5.5 Examples 


5.5.1 Example of a Developable B-Spline Strip Surface 


This example is a windshield of an automobile as depicted in Figure 5-8. Since automobile 
windshields are generally sandwiched between two parallel planes, we can use the technique 
introduced in Section 5.2. The design curve, r4(u), is given in the xy-plane by a cubic 
B-spline curve with control points (0, 2.1, 0), (0.5, 0.9, 0), (1.5, 0.3, 0), (4.5, 0, 0), (7.5, 0.3, 
0), (8:5; 0:9, 0); (9, 2.1, 0) and a knot vector (0; 0, 0, 0, 0:3, 0.5, 0:7 1515 1) 1), Thetis 
end boundary conditions are chosen to be bo = (0.2, 3.5, 1.5) and bg = (8.8, 3.5, 1.5). The 
design curve is first split into four cubic Bézier curves and the control points of rg(u) are 
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R(0,0) 


directnx, rg(u) 





(a) Developable B-Spline Surface (b) Plane Development of Surface 


Figure 5-7: Development of B-Spline Developable Surface 


evaluated. The resulting B-spline curve, rg(u), has internal control points (0.56, 2.64, 1.5), 
(0.96, 2.18, 1.5), (1.84, 1.81, 1.5), (3.27, 1.62, 1.5), (4.5, 1.53, 1.5), (5.73, 1.62, 1.5), (7.16, 
1.81, 1.5), (8.04, 2.18. 1.5) and (8.44, 2.64, 1.5), with knot vector (00000 3 .3 .5 .5 .7 
.71111 1). Note that the resulting surface is not cylindrical. After degree elevating the 
design curve by one and removing the knots from rg(u), we obtain the degree (4-1) B-spline 
developable surface shown in Figure 5-8. The windshield surface is then developed onto the 
plane as shown in Figure 5-8(c). 


5.5.2 Example of a Developable B-Spline Surface with 3D Directrices 


The stack of a ship, depicted in Figure 5-9, was constructed using a B-spline design curve 
with the resulting curve constrained to a plane not parallel to the plane of the design curve. 
Therefore, the method described in Section 5.3 was required. The stack was split along the 
plane of symmetry (the zz plane) from forward to aft. The simple shape allowed definition 
of the design curve with a single B-spline curve of two Bézier segments for each half stack. 
A cubic design curve was selected that defined the base of the stack, with control points (0, 
mem, 10, 0), (15, 10, 0), (35, 3; 0) and (35, 0, 0) and knot vector (0, 0, 0, 0, .5, 1, 1, 1, 
1). The two endpoints were selected to be bo = (6, 0, 31) and bg = (24, 0, 29). The second 
directrix was constrained to lie on a plane, with the third point to define the plane being 
(24, 10, 29). The optimization process was conducted using equations (5.34) and (5.35) to 
reduce the number of variables, and no simple bounds were employed. The resulting surface 
was degree (3,1). The maximum Gaussian curvature determined using equation (5.39) with 
10 Bézier subpatches was 2.19 x 10~’. The surface was then reflected with respect to the 
xz plane, resulting in a symmetric body. One of the surfaces surfaces was developed onto 
a plane, as shown in Figure 5-9(c). 
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(a) Developable Surface with Control Points 


(b) Final Developable Surface 


X 


(c) Surtace Developed onto Plane 


Figure 5-8: Automobile Windshield, a Developable Degree (4-1) B-spline Strip Surface 
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(a) Plan View with Control Points 


(b) Sheer View 
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(c) Half Stack Developed onto Plane 


Figure 5-9: Ship Stack, A Developable Degree (3,1) B-Spline Surface with |K| < 2.19 x 107’ 
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Chapter 6 


Geodesics on Developable Surfaces 


6.1 Introduction 


Geodesics are useful in ship design and shipbuilding for determining the layout of seams and 
butts in a ship hull, as described by Munchmeyer and Haw [33]. This chapter shows that a 
two point boundary value problem (BVP) for computing geodesics on a developable surface 
can always be reduced to an initial value problem (IVP), with a corresponding savings in 
time for the solution process. Most problems that arise in applications of geodesics are 
BVPs, which are much more difficult to solve than IVPs. 


6.2 Formulation 


Recall from Section 2.3 that the governing equations of a geodesic are given by a set of 
coupled second order ordinary differential equations (2.29) and (2.30) 


2 d 2 2 
as +Th (=) oh eee eae (=) = 


ds 12 ds ds ds 
d?y du\? du dv dv\? 
ds? aoa ds see 12 ds ds pene ds Y 


where u and v are related by the first fundamental form ds* = Edu? + 2Fdudv + Gdv?. 
These two second order differential equations can be rewritten as a system of four first order 
differential equations 


a a (6.1) 
= 4 (6.2) 
? = —-[}ip* — 2P}2pq — T9q7 (6.3) 
a = -Dq\p* — Vyopq — Toq°. (6.4) 


We can solve the system as an initial-value problem (IVP), where all four boundary 
conditions are given at one point, or as a boundary-value problem (BVP), where the four 
boundary conditions are specified at two distinct points. Most of the problems that arise 
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in applications of geodesics are BVPs, which are much more difficult to solve than IVPs. 

The solution of an IVP is unique; however for a BVP it is possible that the differential 
equations will have many solutions or even no solution [25]. General methods for the 
solutions of two-point BVPs can be found in [12, 25]. There are two commonly used 
approaches to the numerical solution of a BVP, namely the shooting method and the finite 
difference method (relaxation method). The shooting method was originally used to adjust 
the settings of artillery equipment to hit a target. To implement the method, the unknown 
boundary conditions are guessed at the initial point and an IVP method is implemented. 
Based on the resulting values at the end point, the initial guess is adjusted until the “target” 
is “hit”. The second approach is based on a finite difference approximation to & where 
y = (u,v,p,q) on a mesh of points in the interval [A,B]. This method starts with an 
initial guess and improves the solution iteratively. Maekawa [30] solves the geodesic BVP 
on free-form parametric surfaces by the finite difference method. In general, the shooting 
method is very sensitive to the initial guess at point A, and is much less stable than the 
finite difference method [30]. 

For a B-spline developable surface patch there is a single unique solution to the system. 
Since a geodesic on a developable maps to a straight line on a developed plane, there is 
only one solution to the system on a developable surface. Here we exclude periodic surfaces 
such as cylinders where there can be more than one solution, as shown in Figure 6-1. 





Figure 6-1: Geodesics on a Closed Surface 


In this section it is shown that all two point BVPs for solving geodesics on developables 
can be reduced to IVPs using the following properties described in Section 4.2: 


e A developable surface can be mapped isometrically onto a plane. 
e Isometric surfaces have the same Gaussian curvature at corresponding points. 


e Corresponding curves on these surfaces have the same geodesic curvature at corre- 
sponding points. 


e Every isometric mapping is conformal; i.e. the angle of intersection of every arbi- 
trary pair of intersecting arcs on a developable surface is the same as that of the 
corresponding inverse image in the plane at the corresponding points. 
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e A geodesic on a developable surface maps to a straight line in the plane. 


The basic procedure is to map the two desired points on the developable surface to a plane, 
draw the straight line between them and determine the angle between the generator and 
the geodesic line at one of the end points. The angle can be used to determine du/ds and 
du/ds. Thus, all the information required for an IVP is available. 








C 
(a) Surface Developed onto Plane (b) uv Parametric Space 
B 
g(S) 
se 
A 6 a 
R(0,0) C z 
(c) Geodesic on 3D Surface (d) Final 3D Surface with Geodesic 


Figure 6-2: Geodesic on a Degree (3,1) Developable Surface 


Given two points A and B on the developable surface R(u, v) as shown in Figure 6-2(c), 
the corresponding points (X 4, Y4) and (Xg, Yg) in the developed plane are required. To find 
the points in the plane, (Xo, Yo) is set as the (0, 0) point in the plane corresponding to R(0, 0) 
on the surface. Using the method presented in Section 5.4, the directrix corresponding to 
v = 0 is developed into the plane to determine the point C = (u,4,0), shown in Figure 6- 
2(a). The angle between the directrix and the generator at (w4,0) is the same in both 
representations and can be found by 


Ru (ua, 0) CA 


cos @ = ——————__ - 


FRu(wa,0)] ie) 


where CA is an isoparametric line R(u,, v) which is a straight line on the surface. Therefore 


it is a geodesic and will be developed into the plane as a straight line. The distance ICA 
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is given by 
ICA) = (ca — 20)? + (ya = ye)? + (24 — 2¢)?. 


The point A on the plane is found using C, ICA and @. The point B on the plane is found 
by following the same procedure, and the points are connected as shown in Figure 6- 2( a). 


The angle w between CA and Ale is given by 


This angle w is preserved between the isoparametric line r(w4,v) and the geodesic curve 
g(s) on the developable surface at point A. Thus we have 


R, -g'(s) 
|R»||g’(s)| 


where the tangent vector to the geodesic is given by 


COS W = 


(0) = Rt +R,” 
g (s) =Rus + Ru. (6.5) 


The angles w and @ are shown in Figures 6-2(a) and (c). Multiplying (6.5) by R, yields 
; du dv ; 
Ry, -g(s) = Ru Ry + Ry Ry = cos(w) (Rug (s)| 


which (since |g’(s)| = 1) can be reduced to 


du dv 
F— + G—. = cos(w) VG 


where F = R,-R, and G=R, -R, (coefficients of the first fundamental form). Thus, 





dv _ cos(w) F (=) (6.6) 
ds VG Ne Gas) 
From the first fundamental form, 
du\? du dv dv\? 
/ / a ee mae = 1, : 
g’(s) g'(s) = E (=) + 2F +a(2) 1 (6.7) 


Plugging (6.6) into (6.7) and solving for du/ds yields 


du is sin?(w)G 


de men BG F2 


and thus (6.6) reduces to 





Hence, we have all the initial conditions required to solve the IVP (equations (6.1) to (6.4)) 
for a geodesic. The solution to the IVP yields the uv parametric values for the geodesic that 
are graphed in Figure 6-2(b). The corresponding three-dimensional coordinate values are 
shown in Figures 6-2(c) and (d). The developable degree (3,1) surface shown in Figure 6-2 
has the control points and knot vectors shown in Table 6.1. The geodesic runs from 


(Uas04 |= (0.1, 0.3) fo (ue, Up) 10-9, 0:3). 
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Table 6.1: Control Points for a Degree (3,1) Developable Surface 






6.3 Examples 


In this section we conduct numerical experiments to compare the performance of the IVP to 
the finite difference based BVP on developable surfaces. For a fair comparison, we determine 
the minimum number of steps in the BVP to reliably meet the given tolerance « = 107° 
for Newton’s method, see [31]. We then find the minimum number of steps for the IVP 
to achieve the same order of accuracy at the endpoint B. The examples were run on a 180 
MHz Silicon Graphics workstation. 

The example shown in Figure 6-2 runs using an IVP solution with 66 steps in 0.25 
seconds, as compared to 0.60 seconds for a BVP solution with 135 steps and 7 iterations. 
This is over twice as fast, with a very simple surface and a geodesic that is very close to a 
straight line in uv space. 

A more complex degree (4,1) developable surface with five Bézier patches also runs over 
twice as fast: 0.88 seconds for the IVP solution with 194 steps instead of 2.28 seconds for 
the BVP solution with 294 steps and 11 iterations. The control points and knot vectors for 
this example are displayed in Table 6.2. The geodesic runs from (u,,v4) = (0.1, 0.9) to 
(up,vg) = (0.9, 0.1). The surface and its geodesic are displayed in Figure 6-3. 

The final example is a handkerchief-like surface. This degree (4,1) developable surface 
with four Bézier patches runs in 0.50 seconds for the IVP and 2.66 seconds for the BVP, 
which is over 5 times faster. The IVP required 176 steps, while the BVP required 320 
steps and 12 iterations. The control points and knot vectors for this example are shown in 
Table 6.3. The geodesic runs from (u,,v4) = (0,0) to (wg, vg) = (1,1). The surface and 
its geodesic are displayed in Figure 6-4. 
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b) uv Parametric Space 
(a) Surface Developed onto Plane (b) P 


(c) Geodesic on 3D Surface 


Figure 6-3: Geodesic on a Degree (4,1) Developable Surface 


(0.00, 3.00, -5.00 
(2.42,-1.84, -5.00 
(3.88,-1.53, -5.00 
(5.85, 1.34, -5.00 
( 
( 
( 


) 
) 
) 
) 
(oA) 
8.21,-1.97, -5.00) 
9.09,-2.20, -5.00) 
(9.86,-2.03, -5.00) 
(1.04,-1.53, -5.00) 
(1.09,-6.81, -5.00) 
(1.13, 4.04, -5.00) 
(1.18, 1.05, -5.00) 
(1.20,-1.11, -5.00) 





Table 6.2: Control Points for a Degree (4,1) Developable Surface 
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b) uv Parametric Space 
(a) Surface Developed onto Plane ) 


(c) Geodesic on 3D Surface 


Figure 6-4: Geodesic on a Degree (4,1) Handkerchief-like Developable Surface 


(0.50, 3.00, 0.00) (0.50, 3.00, ) 
(0.88, 1.50, 0.00) (0.88, 1.50, -5.00) 
(1.25, 2.00, 0.00) (1.25, 2.00, -5.00) 
(1.96, 4.17, 0.00) (1.96, 4.17, -5.00) 
(2.50, 3.00, 0.00) (2.50, 3.00, -5.00) 
(3.00, 1.67, 0.00) (3.00, 1.67, -5.00) 
(3.50, 3.00, 0.00) (3.50, 3.00, -5.00) 


(4.04, 4.17, 0.00) (4.04, 4.17, -5.00) 
(4.75, 2.00, 0.00) (4.75, 2.00, -5.00) 
(5.13, 1.50, 0.00) (5.13, 1.50, -5.00) 
(5.50, 3.00, 0.00) o | (5.50, 3.00, -5.00) 
0, 0, 0, 0, 0, 
25, .25, .5, .5, .75, .75, 
lesen 





Table 6.3: Control Points for a Degree (4,1) Handkerchief-like Developable Surface 
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Chapter 7 


Lines of Curvature on Developable 
Surfaces 


7.1 Introduction 


This chapter discusses the behavior of the lines of curvature on a developable surface near 
the inflection line. A developable surface does not possess isolated generic flat points but 
rather may contain a line of non-generic flat points along the generator. This property can 
be used to detect an inflection line where the (only) non-zero principal curvature whose 
principal direction is perpendicular to the generator changes its sign. Munchmeyer [33] 
indicates that a developable surface can be shaped purely by rolling and should be fed into 
the roller so that the direction of the generator is parallel to the rolls. However, when the 
non-zero principal curvature changes sign, the plate can no longer be fed into the roller 
in the same direction, since the bending direction changes at the inflection line. Locating 
the inflection line prior to fabrication allows the manufacturer to either cut the plate along 
the inflection line or reverse the direction the plate is fed through the rollers at that line. 
Surface inflection was studied by Hoitsma [21], who showed that a surface has an inflection 
at a point P if and only if its mean curvature changes sign in the neighborhood of P. This 
chapter further extends this result. 


7.2 Inflection Line 


This section establishes the definition of an inflection line on a developable surface. First, 
some additional differential geometry properties are required. 

Recall that throughout this thesis, it is assumed that the u = const isoparametric line 
corresponds to the generator of the developable surface or, in other words, the straight line 
ruling is in the v direction. With this assumption, R,, = 0; consequently, from (2.13), 
N=R,,-N=0. From (2.23), since Gaussian curvature of a developable surface is zero, 


LN — M? —M? 


oS Se eS = 
EG-F* EG-F? 


0 (7.1) 


and we have M =0. Therefore, mean curvature as expressed by equation (2.24) reduces to 


_2FM-EN-GL_— -GL 


=~ Q(EG— PF) ~ (6G F%)’ oo 
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Recall from (4.22) that the nonzero principal curvature is given by «* = 2H. When L = 0, 
H and x«* become zero; otherwise, k* = 2H £ 0. 

Next, we will show that the u = const parametric straight lines become the lines of zero 
curvature. This can be seen from the fact that 


e the u = const parametric straight lines have zero normal curvature, and 
e no other direction has zero normal curvature. 


The second fact comes from Euler’s theorem (2.25), 
Kn = K1 coS*(a) + Ko sin*(a). 


When «2 = 0, Kn = &1 cos*(a@), which becomes zero only when a = SOE SF corresponding 
to the direction of Kg. Similarly, when x; = 0, «, = 0 only when a = 0 or 7, corresponding 
to the direction of K1. 


Theorem 7.2.1 A developable surface does not possess generic isolated flat points | but 
rather may contain a line of non-generic flat points along a generator. 


Proof: From equation (7.2), L vanishes at a flat point (us, vy) where K = H = 0. There- 
fore, from equation (2.14) we have 


L(us, vs) = —Ry(uy, vs) -Nu(uy, vs) = 0. (7.3) 
From equation (7.1), M@ =0 on a developable surface. Therefore, from (2.15) we have 
M(uy, vz) = —Ry (us, vs) -Nu(us, vg) = 0. (7.4) 
Since N is a unit vector, we also have 
N(uy, vs) -Nu(uy, vs) = 0. (7.5) 


If N, is not zero, then from (7.3), (7.4) and (7.5) N, must be perpendicular to R,, R, 
and N. This is impossible because N is perpendicular to both R, and R,, and R, is not 
parallel to R,. Thus, N,(uys, vs) must equal zero. 

For developable surfaces, the unit normal vector N is constant along a generator. There- 
fore, the rate of change of the unit normal vector in the u direction must also be constant 
along a generator. This leads us to the fact that N, is not only zero at (uys,vys) but also 
zero along the u = uy isoparametric line. Therefore, for a given u = uy, equation (7.3) 
becomes 


L(us,v) = —Rauu(uy,v) > N(uy,v) =0 (7.6) 


for 0 < v < 1. Consequently, the entire generator consists of a line of flat points. I 

For a developable surface, the inflection line is a generator which consists of a line of 
flat points and coincides with a change in sign of the nonzero principal curvature. In other 
words, the nonzero principal curvature has a different sign on either side of an inflection 
line. 


I : *y" . : : « 
A developable surface can not possess spherical umbilics since one of the principal curvatures is always 
zero. 
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The inflection line can be detected by finding u = uy such that L(uy,v,) = 0 where vy, 
is an arbitrary constant between 0 and 1. L(u,v,,) = 0 can be written as 


Ri(u, vn) x R,(u, Un) 
Ryu (u, vp) > = 0. Tf) 
v(t 0m) “TRE (tn) % Rest Un) | 
Since we are assuming a regular surface such that |R, x R,y| 4 0, we need only set the 
numerator of equation (7.7) to zero. Thus, 


Reuss aa, vp) R,(u, Un) = 0. 


For a polynomial surface with degree n in the u direction, this results in a univariate 
polynomial equation of degree (3n — 4) in u 


Yrs = ea = eae = Den One le enue = Dn ae = 0. (7.8) 


If the surface is expressed in a piecewise polynomial form such as a B-spline representation, 
equation (7.8) must be applied to each polynomial segment separately. The univariate poly- 
nomial equation can be robustly and efficiently solved by the Interval Projected Polyhedron 
algorithm {24, 39]. Once we find the u value, we check if the nonzero principal curvature 
changes sign. If the sign does not change, the generator is not an inflection line. A simple 
example is given by R(u,v) = (u,v, u*). Since 


(—4u°,0, 1) 


G=1, u = (0,0,12u7) and N= 
Ru ( ) JVi6u® +1 


we have 


P 1202 


(ee 
V16u® +1 


It is apparent that «* becomes zero when u = 0; however «* does not change sign when u 
moves from u <0 to u>0. Therefore, u = 0 is not an inflection line. 

At flat points the principal directions are indeterminate and the orthogonal net of lines 
of curvature may have singular properties. In the following we investigate the pattern of 
the lines of curvature near the line of non-generic flat points. Locally any surface is the 
graph of a differentiable function [8]. Given a point p of a surface S we can choose the 
coordinate axis so that the origin O of the coordinate is p and the z axis is directed along 
the positive normal of S at p; thus, the zy plane coincides with the tangent plane of S at 
p. It follows that a neighborhood of p in S can be represented in the form z = h(z, y) with 
h(0,0) = h,(0,0) = h, (0,0) = 0. 

By taking into consideration that h(0,0) = h,(0,0) = hy(0,0) = 0, the Taylor expansion 
of the z component of the surface becomes 


1 
na, y) = at hea(0, 0) as 2ryhzy(0, 0) =r yah! 0)] 
1 
+3 [2°heee(0; 0) or Baye O, 0) a Sry" hry, (0, 0) aa y” yyy (0, 0)] 
+R(z,y)(\z,yl°) (7.9) 


where R(z, y) is a remainder term with limz-,0,-40 R(z, y) = 0 and |z, y| = Vz? + y?. 
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This local approximation will now be applied to developable surfaces. 


Lemma 7.2.1 A developable surface is, in general, locally a parabolic cylinder and becomes 
a cubic cylinder at the inflection line. 


Proof: Let us consider an orthogonal Cartesian reference frame attached to the surface at 
an arbitrary point (uo, vo). We choose unit vectors 1] x N, | and N as the directions 
of x, y and z axes such that the y axis coincides with the generator u = ug, the z axis 
coincides with the normal vector and the z axis is orthogonal to both axes. Equation (7.9) 
reduces to [32] 


1 1 
A(z, y) = 5t heo(0, 0) + gt Aaee (0; 0) + Soeur: 0)] + R(x, y)(\x, y|°). (7.10) 
Using the inverse function theorem, hzz(0,0), Arrz(0,0) and hzzy(0,0) can be obtained as 
a function of u and v. At a the line of inflection, hzz(0,0) and Azzy(0,0) become zero [32], 
and equation (7.10) reduces to 


3 
h(z,y) = =(N Ria) ( — x + R(z,y)(|z, yl?) (Gielli.) 


where the coefficient of z° is evaluated at (uo, v9). From equation (7.10) it is apparent that 
for small z, the quadratic term dominates. From equation (7.11) it is apparent that at an 
inflection line, the surface is locally a cubic cylinder. I 


Theorem 7.2.2 There is only one line of curvature that passes through each flat point on 
an inflection line, and that line of curvature is orthogonal to the direction of the generator. 


Proof: By Lemma 7.2.1 the developable surface is expressed locally as a cubic cylinder. 
If we rewrite the cubic term of equation (7.11) in terms of polar coordinates by substituting 
x =rcos@ we obtain 


3 Ee : 
h(@) = = (N : Bare) ( —_ cos? 6. (ale 


Since h(@ + 7) = —h(@), h(@) is an antisymmetric function of 6. The roots of 4 = 0 will 
give the angles where local maxima and minima of h(@) may occur around the flat point. 
The equation can be restricted to the range 0 < 8 < 27 without loss of generality. The 
roots are easily computed as 0 = 0, 5, m and aT Only 6 = 0 and @ = =z (which coincide 
with the local z axis) give extrema, since 


Pn(z) @h(¥) 
de2——s—“‘itséé«?. 


and thus @ = 5 and@ = = (which coincide with the local y axis) provide neither a maximum 
nor a minimum. In other words, the inflection line is not a line of curvature. Consequently, 
there is only one line of curvature (minimum curvature where 0 = 0 and maximum curvature 
where @ = 77) that passes through the flat point. I 
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Recall from Section 2.3 that the principal curvatures can be found by solving the two 
simultaneous equations (2.19) and (2.20), 


(L+«E&)du + (M + «F)dv 
(M+«F)du+(N+«G)dv = 0, 


for «. Since M = N = 0, equations (2.19) and (2.20) reduce to 


(L+«F)du+ «Fdv 
KFdu + KGdvu 


0 and (7.13) 
0. (7.14) 


We can trace the lines of curvature which pass through the flat points of an inflection line 
by providing the starting points and shifting outwards in the directions 0 and 7 from the 
flat points or, equivalently, along the positive and negative local z axis. The tracing can 
be carried out’ by integrating the initial value problem for the coupled ordinary differential 
equations (7.13) and (7.14). 


Generic umbilics (including flat points) are stable with respect to small perturbations of 
the function representing the surface, while non-generic ones are unstable [32]. In generic 
cases umbilics are isolated [20]; thus the inflection line, which consists of line of flat points, is 
non-generic and therefore unstable. In the following we give a couple of numerical examples 
that demonstrate the instability of the line of flat points along the inflection line with respect 
to perturbations. 

The example surface is a degree (3-1) integral Bézier patch which is constructed by the 


method described in Section 5.3, see Figure 7-1. The control points are shown in Table 7.1. 
The surface has an inflection line at u = 0.5754, which has been computed by solving the 


0.0000, 0.0000, 0.0000 
1.8000, 3.0000, 0.0000 


(0.5000, 0.0000, 2.0000) 
(1.8950, 2.3250, 2.0000) 
(3.0575,-1.5500, 3.1625) 
(3.6000, 0.0000, 2.0000) 





) 
) 
3.3000,-2.0000, 1.5000) 
4.0000, 0.0000, 0.0000) 


Table 7.1: Control Points for a Degree (3,1) Developable Surface 


degree 5 univariate polynomial equation (7.8). The inflection line is shown in Figure 7-1 
with a dash dotted line. This surface has a net of lines of curvature which is shown in 
Figure 7-2(a). Solid lines represent the lines of maximum principal curvature, while dotted 
lines represent the lines of minimum principal curvature. The inflection line is depicted with 
a dash dotted line. Figure 7-2(b) shows a magnification near the inflection line. We can 
observe that there is only one line of curvature that passes through a flat point orthogonal 
to the inflection line. 


We gradually perturb the control points of the surface and observe the behavior of the 
lines of curvature which pass through the inflection line. The control points are perturbed in 
the following manner. Since the example is a degree (3-1) patch, it has 8 control points. Each 
control point consists of three Cartesian coordinates x, y, z; hence there are 24 components 
to be perturbed. A random number which varies from —1 to 1 is used to determine the 
perturbation of each component. Let us denote the randomly chosen numbers for each 
control point as (ez, er, oF 0 <1i,7 <3. We normalize the vector and add it to each 
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Figure 7-1: Developable surface and its control polyhedron with inflection line. 





inflection line 
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Figure 7-2: (a) Lines of curvature of developable surface with inflection. (b) Magnification 


near inflection line. 
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Figure 7-3: (a) Lines of curvature on perturbed surface ¢ = 0.02. (b) Magnification near 
u=0.97. 


_- PF eM 


—-_ — 
—_—|— 
= 
_ 
_ 





\ 
~ 


Figure 7-4: (a) Lines of curvature on perturbed surface ¢ = 0.08 (b) Magnification near 
u=0.07. 
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control point as follows: 


Yo ez. 
Jd (7.15) 


/ 22 y 2 Zz 2 
Cr; T €;; + &;, 


where ¢ = 0.02. We gradually increase the perturbation by increasing ¢ from 0.02 to 0.08 
in steps of 0.02. | 

Figures 7-3 and 7-4 illustrate the behavior of the lines of curvature when the control 
points are perturbed. We can see from the figures that the entire inflection line, which 
consists of a line of flat points, disappears. Hence there is no singularity in the net of lines 
of curvature when a perturbation is induced. The nonzero principal curvatures on both 
sides of the former inflection line 7 meet at right angles near the former inflection line and 
make a very sharp change in direction (almost a right angle). 


i [e7., e 
P;; = P3; as ¢ J 


*Once the control points are perturbed both principal curvatures may not be nonzero, but here we are 
referring to the nonzero principal curvature before perturbation. 
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Chapter 8 
Engineering Example 


The final example is the design of a small boat using B-spline developable surfaces. The 
boat was designed with a hard chine which separates the bottom and side of the boat and 
demarks a tangent discontinuity in the hull of the boat. As an initial stage of the design, the 
following three curves, which define the desired characteristics of the hull, were determined 
using a model of a planing boat: 


e sheer curve — the line that runs along the deck edge from bow to stern 


e base curve - the line that runs along the keel of the boat (centerline from bow to 
stern) 


e chine curve — the line that runs from the bow to the lower corner of the stern 


Figure 8-1 shows the terminology used in this section. The coordinate system is located 
so that the origin is located on the baseline at the bow of the boat with the z axis pointing 
aft, the y axis pointing to starboard and the z axis pointing up. (Note that the baseline, 
a zero reference line, is different from the base curve defined above.) Due to the symmetry 
of the boat, we design only the starboard side and reflect those values across the plane of 
symmetry to obtain the port side. The cubic B-spline representation of the three curves is 
shown in Table 8.1. 


Transom Chine Sheer 


Side Aft "es Va 
“a 7 
co Bottom Aft ae 
Sg cree — y 


Figure 8-1: Planing Boat Terminology 














Design Using 3D Directrices 


Employing the design philosophy described in Chapter 5, the chine curve was designated 
as the design curve and the sheer and base curves as target curves, i.e. rg(u), to construct 
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[Sheer Curve [Chine Curve [Base Curve 

























Control | (0.00, 0.00, 9.00) | (1.40, 0.00, 5.30) | (1.40, 0.00, 5.30) 
Points | (6.86, 7.10, 8.22) | (10.5, 7.53, 1.93) | (2.26, 0.00, -0.21) 
(21.6, 8.93, 6.25) | (25.7, 7.85, 1.28) | (22.6, 0.00, -0.10) 
(36.9, 8.73, 5.86) | (40.4, 7.46, 1.27) | (36.3, 0.00, -0.10) 
(45.0, 7.65, 6.10) | (44.1, 7.20, 1.70) | (44.1, 0.00, 0.50) 





Knot ae 20, 0, .o, CUO Onto, OOO Ome. 
Vectors tol ae ele lead eel al 


Table 8.1: Target Curve Data 


the side and bottom respectively. Initially, the endpoints of the sheer and base curves were 
used to construct the side and bottom of a single surface each; however, this resulted in a 
flat side and bottom. To remedy this situation, the chine was split into two B-spline curves 
at a parameter value of u = 0.5, and one additional knot was inserted at the midpoint of 
the forward B-spline curve. The resulting curves are shown in Table 8.2. Thus, the side 
and bottom are constructed of two surfaces each, allowing more flexibility in the design, 
but losing the C* continuity along the surface. The surfaces are, however, at least G! 
continuous between the forward and aft portions. Since the design curve was split into two 
curves that maintain C? continuity between them, the tangent at the connecting point is 
continuous. Since both patches are developable, both surfaces share the same tangent plane 
at the connecting generator. Hence, the patches are G! continuous. 


Note that the forward and aft surfaces join along a generator which is a straight line. 
The planar transom (stern) joins the aft surfaces along generators as well. 


The endpoints for the side and bottom surfaces were selected so as to force the resulting 
directrix as close to the target curves as possible while maintaining a maximum Gaussian 
curvature on the order of 10~°. The data for each developable surface on the starboard side 
of the boat is listed in Table 8.2. Measurements are in feet. The port side of the boat has 
the same measurements, with a negative y value. 


Patch Design Curve Knot Endpoints Plane Gaussian 
| contin | ver | "P| peinion ot | Cate 
Forward 1.40, 0.00, 5.30 ae 0.00, 9.00) | (20.0, 0.00, 6.74) | 8.18 x 107 
Forward 10.7, 5.86, 2.56 ek 5 10, 0. 00, OD. 30) 
Bottom | (18.2, 7.09, 1.87 (14.0, 0.00, 0.00) 


(0.00, 0.00, 0.00) 
21.9, 7.38, 1.66 
Aft 21.9, 7.38, 1.66 (20.0, 8.30, 6.74) | (45.0, 0.00, 6.10) | 6.62 x 10— 
Aft 33.0, 7.66, 1.28 14.0, 0.00, 0.00) | (0.00, 0.00, 0.00) | 9.74 x 10- 
Bottom | (40.4, 7.46, 1.27 44.1, 0.00, 0.50) 
44.1, 7.20, 1.70 





Table 8.2: Developable Surface Data, Starboard Side 
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Constrained by Plane 


All the second directrices were constrained so that they lie in a plane defined by the two 
designated endpoints of rg(u) and a third point g provided by the user. This point g is 
termed the Plane Definition Point and is shown in Table 8.2. For the bottom surfaces, 
the constraining plane was the xz-plane which divides the boat into symmetric halves. For 
the side surfaces, the constraining plane was positioned so that a line in the plane and 
perpendicular to the z axis would be parallel to the zy plane. 

Figure 8-2 shows the resulting developable surfaces with the design lines shown as dotted 
lines for reference. Note that the forward portion of the bottom of the boat is a triangular 
degenerate patch. Figure 8-3 is a bottom view of the boat. One side shows the rulings of 
the developable surfaces while the other side is a shaded view. 





Figure 8-2: Boat Composed of Developable Degree (3-1) B-Spline Surfaces 





WZ 
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Figure 8-3: Boat Bottom View 


ol 


Development 


The plane development of each surface of the boat was computed, thus providing cutting 
information for the flat steel. Using the method described in Section 5.4, the planar surfaces 
shown in Figure 8-4 were constructed. All developments were within a tolerance of 107°. 


[a ay eae 


(a) Aft Side (o) Forward Side 


(c) Aft Bottom (d) Forward Bottom 


Figure 8-4: Boat Surfaces Developed onto a Plane 


Geodesics 


At the Electric Boat Corporation shipyard, many developable surfaces are used in the 
construction of submarines, since a large portion of the hull is cylindrical. When bending 
the steel into shape, accuracy is checked with detailed measurements of geodesics drawn on 
the surface. When the surface is flat, straight lines are drawn between the corners of the 
surface, forming an X as shown in Figure 8-5(a). The position of this line on the curved 
surface is easily computed using the method described in Chapter 6. As an example, the 
geodesics on the Forward Side section have been calculated. The uv parameter values are 
shown in Figure 8-5(b) and the geodesics on the developable surface are shown in Figure 8- 
Oe} 


Lines of Curvature 


These surfaces have no inflection lines, but the lines of curvature were calculated as described 
in Chapter 7 and are shown in Figure 8-6. The dashed lines represent lines of minimum 
curvature and the solid lines are lines of maximum curvature. Note that one set is always 
composed of straight lines that correspond to the generators. 
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(a) Surface Developed onto Plane (b) uv Parametric Space 


(c) Geodesic on 3D Surface 


Figure 8-5: Geodesics on the Forward Side Section 
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Soaeeeeeeey 


(b) Forward Side a) Aft Side 
ye 
ERS oe fp 
(d) Forward Bottom (c) Aft Bottom 


Figure 8-6: Lines of Curvature on Boat Surfaces (Solid Lines are Maximum Curvature 
Lines, Dashed Lines are Minimum Curvature Lines) 
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Chapter 9 


Conclusions and Recommendations 


9.1 Conclusions 


In this thesis a novel method to design complex objects in terms of B-spline developable 
surfaces focusing on user friendliness is developed. The user only needs to specify the 
design B-spline curve and the two end rulings, and the developable surface in a B-spline 
representation is automatically generated. Although optimization techniques are relied 
upon, in most cases the computed results have a Gaussian curvature on the order of 107° 
or less unless the design of a very complex surface is attempted. The computational time 
for the design of all the examples in this thesis is on the order of a few tenths of seconds on 
a 100 MHz Silicon Graphics workstation. 

In addition, two differential geometry properties of developable surfaces which have 
practical engineering applications were investigated. The first property was the definition 
of inflection lines on developable surfaces, which are critical in the manufacture of such 
surfaces. We proved that an inflection line on a developable surface occurs on a generator 
which consists of entirely flat points, and derived a single variable nonlinear polynomial 
equation to locate it. It was also found that for each flat point on an inflection line of a 
developable surface, there is only one line of curvature that passes through it, and that line 
of curvature is orthogonal to the direction of the generator. 

The second differential geometry property investigated was the location of geodesics 
between two points on a developable surface. It was shown that the problem can always be 
converted from a boundary value problem to an initial value problem with a corresponding 
savings in time for the solution. 

Finally, these techniques were applied to a boat model, a complex object that was 
represented entirely with B-spline approximate developable surfaces with cubic directrices. 
The Gaussian curvature was on the order of 107°. Geodesics and lines of curvature were 
calculated on these surfaces as well. 


9.2 Future Work 


The ultimate goal of this research is to be a useful tool in the design and manufacturing of 
ship hulls. Therefore, future research should take this work in that direction. Some specific 
areas of possible research are described below. 

Approximation. The approximation of free form surfaces by developable surfaces would 
be extremely helpful in converting existing ship designs to ones including more developable 
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surfaces. Some preliminary research has been completed by Elber and Cohen [9] into repre- 
senting surfaces as ruled surfaces, which could possibly be extended to include developable 
surfaces. In addition, research on approximation using developable surfaces was conducted 
by Hoschek and Pottmann [23] using the representation of developable surfaces as the en- 
velope of a family of planes. Alternatively, Cho [6] discussed the approximate development 
of non-developable doubly curved surfaces onto a plane, which may also be useful in the 
above approximation problem. 

Manufacturing. Once the surface is developed into a plane and the steel is cut to the 
‘correct size, it must be bent into the proper shape. Additional research should investigate 
the proper amount of curvature to be induced by rolling or pressing and instructions for 
the shop floor should be output. Although the lines of curvature are instrumental in this 
bending, mechanical considerations such as springback must be included in the design as 
well. Much research has been done in this area by Hardt et al [19, 18, 17]. The ultimate 
goal is for the engineer to design the ship hull and output information that is directly useful 
to the workers on the shop floor. To this end, the operating characteristics of bending 
equipment should be studied and properly accounted for. 

Geometry. Several areas of interest remain for investigation within the differential ge- 
ometry discussed in this thesis. First, a suitable non-dimensional Gaussian curvature ex- 
pression should be developed to better judge the developability of surfaces designed using 
the minimization process described in Section 5.3. 

Although strip surfaces of any degree can be handled by the computer programs devel- 
oped during this research, the programs for the surfaces with arbitrary 3D directrices can 
handle only B-spline developable surfaces with cubic directrices. It may be beneficial to 
investigate higher degree surfaces which may include more flexibility and would be able to 
more closely approximate the desired output surface. 

In addition, the surfaces in this thesis are all integral B-spline surfaces. Extension 
to include rational B-spline developable surfaces such as those investigated by Pottmann 
and Farin [37] and Lang and Réschel [29] would allow the representation of more complex 
surfaces. 

One difficulty in designing surfaces using this method is that, for surfaces with 3d 
directrices, there is no method to determine in advance if the surface will include the edge 
of regression. Therefore, regular surfaces are achieved through trial and error, and hence 
regularity conditions should be investigated more fully. 

Design. Much research has been conducted in design for producibility in ship hulls. For 
example, Lamb [27] shows a method to construct a bulbous bow from regular surfaces and 
discusses the simplification of bow and stern sections to include more developable surfaces. 
This should be expanded to investigate the design of an entire ship hull from developable 
surfaces. Once a ship hull is designed, it should be hydrodynamically tested to compare the 
performance to a similar hull that is not constructed of developable surfaces. The relative 
magnitude of the savings due to producibility and the losses due to performance should be 
compared. 
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